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VER A PERIOD of several years, 
a number of transformers and in- 
duction regulators supplying energy to 
aluminum potlines and similar installa- 
tions have been damaged by insulation 
failures; the failures were apparently 
due to abnormal overvoltages. The 
transformers used for these applications 
usuaily step down from 13.8 kv to a 
secondary voltage in the range 600 to 900 
volts depending upon the installation. A 
basic impulse level of at least 40 kv is 
normal for the secondary winding; in- 
sulation failures occurring on these wind- 
ings are an indication of the severity of 
the voltage transients. An excellent 
photograph showing typical coil damage 
appears in a paper by Pope, et al.? 

Naturally, these overvoltages have 
engendered considerable speculation as 
to their origin and, although the prin- 
‘cipal object of the work to be described 
here is not to determine the precise 
cause of overvoltages, the authors’ specu- 
lations are added, as a help in understand- 
ing how the work was planned. Later, 
we will look at these speculations again 
in the light of the results obtained. 

It is clear that voltage transients must 
be generated within the system or else 
enter as surges from the high-voltage line. 
Many indications, at least in some 
instances, point to the unlikelihood of 
the latter possibility. Bourque and 
Coxe have described, in unpublished 
work, an installation in which the trans- 
formers were supplied on their high- 
voltage side, through considerable lengths 
of underground cable. This certainly 
minimizes the chance of damaging surges 
approaching by that route. Also, pot- 
line installations frequently comprise a 
number of parallel transformers which 
share the energy of impinging surges, and 
it is difficult to see how transient voltages 
of the order experienced could be trans- 
ferred to the low-voltage winding, except 
perhaps as an initial spike by electrostatic 
coupling, without flashover occurring on 
the high-voltage side, For these reasons 
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‘element or another. 


our attention was focussed mainly on the 
possibilities of surges being developed 


within the system. 


When one looks at a network, such as 
a potline-rectifier-transformer system, it 
would seem that the only way in which 
abnormally high voltages could be pro- 
duced internally would be by some form of 
current suppression. Energy stored in the 
magnetic field of the circuit would thereby 
be released and transferred to the elec- 
tric field, where it would enhance the 
voltage across the system capacitance. 
There appear to be four active elements 
in the network that might bring this 
about: the anode and cathode circuit 
breakers, the rectifiers, and the potline 
itself. Several investigators? * have ruled 
out the circuit breakers since significant 
current chopping is not a characteristic 
of the breakers used. The authors con- 
cur in this opinion; if it is true, it seems 
likely that the disturbances originate in 
the rectifiers or the potline by some 
mechanism akin to chopping. Several 
such possible mechanisms have been pro- 
posed by Hull and Elder* and more re- 
cently by Steiner and Strecker,’ and under 
certain conditions surging in rectifiers 
has been observed in controlled test. 
Potline cells are also known to behave 
rather erratically from time to time. 
However, in our investigation we were 
not concerned with mechanisms or, in- 
deed, with fixing responsibilities on one 
Rather, this study 
starts with the assumption that abrupt 
changes of current could occur, and 
proceeds to analyze the consequences 
thereof. 

For a better appreciation of the results 
of this investigation a brief résumé of 
the behavior of a circuit when the current 
in it suddenly changes, would seem to be 
in order. In the limit when the current 
changes instantaneously, the voltage 
generated is given by the product of the 
instantaneous change of current and the 
surge impedance of the circuit. It is not 
the product Ldi/dt, which in this instance 


would be infinite, for this ignores the 
capacitance of the component or circuit 
involved. Most practical circuits com- 
prise a number of elements in series/ 
parallel combinations, but on analysis 
they can usually be grouped into several 
series units each with its own surge im- 
pedance. 

When stich a circuit is excited by a step 
function of current, or by what for con- 
venience we will call a chop of current, 
each unit will respond in the manner 
described above, i.e., it will be set into 
oscillation at a frequency peculiar to 
itself and, in the absence of damping, the 
peak voltage generated will be 7Zo, 
where Z, is its surge impedance. The 
total response of such a circuit will be the 
summation of the individual responses of 
these units, ic., a multifrequency dis- 
turbance. 

In practice it frequently happens that 
one part of the circuit has a much higher 
surge impedance than the rest so that 
voltage tends to pile up across this part. 
Again, in all practical circuits the sup- 
pression of current always takes a finite 
time although in some cases this may 
be very short indeed. Because of this 
the overvoltage is not as great as the 
theoretical limit. The discrepancy will 
depend upon the rate of decline of cur- 
rent and the character of the circuit it- 
self. Specifically, it has been shown by 
Lee® that the overvoltage will be signifi- 
cant if the current collapses in a time less 
than half a cycle of natural frequency of 
the circuit involved. Declining as a 
linear ramp in this time, the current would 
generate an overvoltage of approximately 
60% of the theoretical peak; for longer 
times, the overvoltage is progressively 
less. This point has a very important 
bearing on the results to be described be- 
low. With this background we can now 
outline the aims and scope of the in- 
vestigation. 

If the voltage transients are, indeed, 
caused by the sudden change of current 
at either a rectifier or in the potline, 
then by creating or simulating a chop 
at these locations, it should be possible 
to investigate the magnitude and dis- 
tribution of the voltages that they give 
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rise to. This is tantamount to measur- 
ing dynamically the surge impedances of 
the constituent parts of the circuit for 
various possible circuit connections, and 
thereby serves two purposes: (1) It per- 
mits one to determine the minimum 
magnitude of the current suppression 
that must occur to generate any specified 
voltage at any specified point. (2) If in 
any future investigation, with a surge 
recorder, for example, a certain pattern 
of transient voltages is disclosed, it 
should be possible by comparing the 
pattern with the test results to pinpoint 
the origin of the disturbance. It was 
therefore not the aim of the investigation 
to determine the cause of overvoltages, 
but rather to say that if current sup- 
pression were the cause, it must occur 
in a particular locality with a particular 
severity. 

The opportunity to carry out such a 
test program came in June 1958, when the 
Reynolds Metals Company was about to 
commission a new aluminum potline at 
their Lister Hill, Ala., plant. The re- 
sults obtained apply in detail to that 
particular system only, but in a broader 
way they provide information on many 
similar circuits. Moreover, the test 
methods described below can be applied 
equally well elsewhere. 


Circuit Arrangements and Test 
Procedures 


A schematic diagram of a typical recti- 
fier skid and associated transformer and 
breaker equipment is shown in Fig. 1, 
where it will be noted that the trans- 
former has two secondary windings. 
In fact, it had four such windings and two 
primaries on a cotmmon core ; the others 
served an identical skid adjacent to the 
one shown. In all, there were 7 such 
transformers and 14 rectifier skids supply- 
ing the potline. 
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Fig. 1 (left). Circuit 
diagram of transformer, 
switchgear, and rectifiers 


SAWTOOTH 
GENERATOR 


Fig. 2 (right). Test cir- 
cuit for current injection 
technique 


Two methods of testing were adopted. 
In the first, a direct current was fed into 
the system at the supposed point of sup- 
pression, e.g., at the anode and cathode 
of one of the rectifiers. This current, 
approximately 0.4 amp (ampere) was 
then interrupted suddenly by means of a 
vacuum switch. In the second method, 
square pulses of current were repeatedly 
injected into the circuit at the rectifier. 
The front of each square pulse represents 
a negative chop, as it were, while the back 
represents a true chop. The voltages 
generated by these sudden changes of 
current were then investigated oscillo- 
graphically at all locations of interest. 
The circuit for the current injection 
method is shown in Fig. 2. It has the 
advantage that, because of the recurrent 
nature of the stimulus, the oscilloscope 
Sweep can be synchronized with the surge 
generator, so that the circuit response is 
seen as an apparently steady picture on 
the screen. With the d-c interruption 
method the oscilloscope had to be trig- 
gered at the initiation of the transient, 
which then appeared as a single sweep on 
the screen. ; 

In both methods the systems were dead 
at the time of the tests. It was recog- 
nized that surges might occur under 
many different circuit conditions, viz., 
normal operation with or without com- 
mutation, backfire, misfire, etc.; these 
different conditions had to be simulated, 
which was a fairly simple process. Those 
rectifiers which were assumed to be con- 
ducting at the time the surge occurred 
were simply short-circuited out with clip 
leads; the remainder were left open, 
This reveals a further advantage of the 
recurrent surge method of testing: changes 
from one condition to another could be 
made quickly and simply while the sys- 
tem was being excited; the consequences 
o* any change could be observed at once 
on the oscilloscope screen, 
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The benefits of repetitive surge techi 
nique has been very aptly summarizee 
by Wilkinson:’ a 


> | 
...a recurrent surge oscillograph is strik 
ingly effective in its power of conveying ing 
formation to the user. This is so because on 
the ease with which questions about the 
effect of any particular detail in a network 
under review are put to the test and imp 
mediately produce a clear fixed image on 
the oscillograph by way of reply. If the 
recurrent surge oscillograph is like a tutos 
with whom problems can be discussed, tha 
single-record oscillograph is to be comparee 
with a system of correspondence by field 
post card. | 


When the _ recurrent method ii 
checked against that using the vacuunt 
switch, the agreement was found to bo 
remarkably good. Because of this agrees 
ment and the advantages of the current 
injection method, this method was use¢ 
almost exclusively during the investi gas 
tion. For the majority of the tests tha| 
13.8-kv bus was short-circuited ana 
grounded in order to simulate the norm a2 
energized condition. The source is stiffi 
so, for all practical purposes, we cami 
neglect its impedance. A few tests were 
conducted with the transformer primaries 
open-circuited, a condition which occurs 
when a potline is being de-energized by 
first interrupting the a-c supply. If ei 
surge should occur in the short intervai 
between dropping out the 13.8-ky circuti 
breakers and the subsequent opening 0) 
the anode and cathode breakers, then ii 
would be simulated by the conditiona 
just described. ] 
So far, nothing has been said about tha 
potline itself. It could not, of coursed 
be introduced directly into the proble mi 
Before the test program was started, the 
problem of simulating the potline waa 
given considerable attention. The pott 
line is physically extensive, in some way7 
resembling a transmission line loopet( 
back on itself; it is isolated from grounc| 
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ig. 4 Response at points 7 and 22 to 20 ma 
(milliamperes) injected at tank 2 


nm both polarities; but it has, one would 
ppose, considerable capacitance to 
It was felt that as a first 
pproximation, the simplest analog would 
be a parallel L-C (inductance-capacitance) 
ircuit. By temporarily disconnecting 
me of the operating potlines and shock- 
ng it, it was hoped to establish some 
eliable value for the Z and C compo- 
1ents. 

When the test program preliminaries 
vere under way, it was both interesting 
nd gratifying to find that the transients 
being measured were relatively insensitive 
o the potline analog. For example, put- 
Hing a 0.5-uf (microfarad) capacitor across 
l-c bus had only secondary effects on the 
magnitude and waveform of the voltages 
renerated around the rectifier skid and at 
the transformer, when simulated current 
thops were applied. The reason for this 
will become apparent when some of the 
Fesults are analyzed. 


There were a number of rather unusual 
features about the system under test that 
presented their peculiar problems from 
2 measurement point of view. Perhaps 
he most important of these was the fact 
hat the entire system was ungrounded, 
though there was capacitance coupling 


Taking measurements on such a system 
Wwas rather like standing on a spring 
mattress and trying to measure the 
Histance to the floor. This situation 
balled for the isolation from ground of 
both the surge-generating equipment and 
the oscilloscope used for measurement. 
The equipment that generated the current 
bulses was fed through an isolation trans- 
ormer which had _ significantly less 
rapacitance between its windings than 
had its own power supply transformer. 
By increasing the impedance of this path, 
e errors that would be introduced by 
Mound current circulating through the 
system capacitance to ground and back 
o the generator through its ground capac- 
tance would be minim zed. 
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The measurement end was tackled by 
measuring the transient voltage differ- 
entially, i.e., a double-ended or differ- 
ential amplifier was introduced between 
the point of measurement and the oscillo- 
scope plates. It was apparent from 
oscillograms taken, that the system 
response was quite intimately dependent 
on the stray capacitance of the system. 
To have used a single-ended amplifier 
input to the oscilloscope would have been 
to connect the capacitance of the oscillo- 
scope (i.e., the capacitance across its power 
supply to the 60-cycle system) to one of 
the measuring points. This would in- 
evitably have introduced an error. 

The oscillogram, Fig. 3, is typical of 
the kind of response that was observed. 
The point of interest here is the very fast 
rate of rise of the initial spike of voltage: 
it reaches its crest in approximately 0.3 
usec (usecond). It will be recalled that 
at the start of the paper a special point 
was made of the fact that the surge gen- 
erated by a current chop depends woon 
the rate at which the current declines; 
the more rapid the rate, the greater the 
voltage, being a maximum when the chop 
is vertical. To approach this maximum 
the time for the current to come to zero 
should be short compared with the period 
of the highest frequency excited. 

In view of the obviously very high fre- 
quency that appears in Fig. 3, it was 
necessary to take a careful look at the 
front of the injected current pulse to make 
quite certain that it was fast enough to 
develop a voltage approaching the maxi- 
mum. Fig. 4 is a record of the current 
pulse wavefront; note that the current 
reaches approximately 80% of the 
maximum value and perhaps 90% of 
the final value on 0.2 usec. It is there- 
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Fig. 4. Profile of injected current 


fore safe to assert that the response shown 
in Fig. 3 is not less than 85% of the maxi- 
mum available from a vertical chop. 

The rise time of the initial transient in 
Fig. 3 raised the question as to whether 
the response of the amplifier was adequate. 
This point was checked with some care 
using different amplifiers... Widening the 
bandwidth did not change the appearance 
of the fast transients. 

Several other smaller difficulties had to 
be disposed of before the test procedure 
ran smoothly. The physical distance 
involved caused some inconvenience; 
in many instances the point of measure- 
ment was quite remote from the point 
where the stimulus was applied to the 
circuit. For example, on one occasion, 
measurements were made at the bushings 
on top of the transformer, situated out- 
side in the yard, and the pulsing equip- 
ment was located inside the rectifier 
room close to one of the skids. 

Synchronizing becomes difficult under 
such conditions, the firing circuits of 
neighboring energized potlines put out a 
constant procession of pulses giving a high 
background noise level, which calls for 
careful screening of trigger circuit leads 
to the oscilloscope. It is usual to have 
the pulse generator initiate the oscillo- 
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Fig. 5. A—Circuit diagram for an arc back 
with one Y conducting. B—Equivalent circuit 
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scope sweep somewhat ahead of the pulse 
but in some of the applications where 
greater physical distances were involved, 
it was found more convenient to reverse 
the process, i.e., the pulse generator 
was triggered by a signal emanating 
from the sweep generator of the oscillo- 
scope. 

A minor irritation was occasioned by 
interference with the electron beam 
in the oscilloscope caused by the strong 
d-c magnetic fields set up by the current 
in a neighboring d-c bus. This could be 
reduced somewhat by aligning the beam 
axis as closely as possible with the 
direction of the magnetic field. It was 
ultimately reduced. to negligible propor- 
tions by putting the oscilloscope in a 
welded steel box made from a 1/4-inch 
boiler plate. 


Analysis of Results 


As explained in the section ‘Current 
Arrangements and Test Procedures,” the 
procedure in test was to connect jumpers 
across all tanks assumed to be conduct- 
ing at the instant under consideration, 
then to inject the simulated chop at the 
point of interest. Different circuit ar- 
rangements were tested in turn and for 
each the response was explored, starting 
close to the point of application of the 
disturbance and working away to more 
remote points to where ultimately the 
effects were insignificant. 

In analyzing the results below, this 
same procedure will be followed. Which 
components were contributing most to the 
response became apparent as the program 
progressed. Once this had been definitely 
established, it was possible to reduce 
appreciably the number of tests in later 
series. To illustrate, one series of tests 
will be discussed in detail here, and im- 
portant information gained from other 
series will be discussed briefly. The 
particular condition to be treated is that 
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in which a surge follows an arc-back with 
one Y only conducting. No special 
significance attaches to this condition, but 
it serves well to illustrate the behavior of 
the system, and demonstrates the method. 

It is assumed that the current sup- 
pression occurs in the rectifier that has 
arced back. A rough circuit diagram of 
the condition is shown in Fig. 5(A), 
where the manner in which the injected 
current divides is indicated. Fig. 5(B) 
is an approximate equivalent circuit. 
Points in the circuit are identified by 
the same numbers as those used in Fig. 1; 
the physical location of these points can 
be understood by the layout as shown 
in Fig. 6, 

Consider first the trace observed be- 
tween points 7 and 22, shown in Fig. 3. 
Since 7 and 22 are fairly close electrically 
to the points of injection, this oscillo- 
gram represents the response of most of 
the system. We notice several features: 
in the first place it is multifrequency. 


Second, it reaches its crest in approxi- 


mately 0.3 usec; therefore one of its 
important component frequencies is quite 
high. Third, the peak magnitude is of 
the order of 60 volts/amp of chopped or 
injected current. One might say that this 
corresponds to an effective surge im- 
pedance of 60 ohms; we call this ‘‘effec- 
tive’ because, as has been pointed out, 
the response is really a composite one, the 
summation of several separate responses; 
also, it inevitably includes damping. In 
the next few oscillograms, it will be possi- 
ble to investigate as to what circuit com- 
ponents are giving rise to the different 
constituents, and to estimate their rela- 
tive importance. 

A quick check across points 21 and 22 
indicated that the dividing reactor located 
between the tanks contributed very little. 
On the other hand, it was found that there 
was a remarkable change at location 9B to 
15B, shown in Fig. 7(A). The high- 
frequency component is quite drastically 
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A—At points 9B and 15B 
B—At points 7 and 21 with short circa 
between points 9B and 15B 
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reduced but the lower frequency remains 
virtually unchanged. This is very ini 
formative; it means that the high fre 
quency is generated, for the most partt 
by the comparatively short lengths of bus! 
between the rectifier skid and the base? 
ment wall and the negative disconnec? 
(the B in the designations 9B and 15] 
indicates that the measuring points are 
in the basement; see Fig. 6). 

We conclude that the lower frequency 
is the product of equipment beyond the 
basement wall, namely, the interphase 
transformer and the power transformer 
outside in the yard. This was furthes 
confirmed and pinpointed by subsequeni 
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| ig. 9. Response to 20 ma injected at tank 2 


WA—At points 7 and 8 
B—At points 7 and 21 with 13.8-kv winding 
of transformer open-circuited 
—At points 7 and 8 injected at these two 
points 


Points 9B and 15B were short- 
ircuited, and the response observed across 
points 7 and 21; the effect is shown in 
Fig. 7(B). The low frequency has dis- 
appeared but the high frequency re- 
mains. 


(points 2 and 44) and the response was 
fmeasured at the same points; see Figs. 
(A) and 8(B) with different time scales. 
Wow the lower-frequency component is 
Pvident, together with some high fre- 
Quency; this arises from the length of 
Hus from the transformer to the anode 
breaker, which could not be readily dis- 
bonnected. Fig. 8(B) reveals a condition 
which could not be appreciated from the 
Dreceding oscillograms because of their 
time scales: the presence ofa yet lower 
Krequency component of significant pro- 
jortions. It was thought that this was 
rontributed by the interphase trans- 
former, an assumption confirmed by ob- 
aining the transient across two phases 

points 2 and 6). This component was 
ound to have disappeared. 
| A very significant fact was revealed 
by taking measurements at points 7 and 
8. The voltage generated across the other 
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13 skids of the potline was thereby 
observable. This voltage is shown in 
Fig. 9(A) and is seen to be quite small rela- 
tive to the voltage on the skid and trans- 
former where the surge originates. This 
is not surprising because the other 13 
skids are connected in parallel through 
the d-c bus and must therefore present a 
low surge impedance. It means that 
when a current suppression occurs, its 
effects are localized to the transformer 
and rectifier equipment immediately 
involved. 


To recapitulate the results thus far, 
it may be said that a surge in a rectifier 
under the conditions prescribed excites 
three principal components. The highest 
in frequency, approximately 1 mc (mega- 
cycle), and amplitude, is caused by the bus 
run from the rectifier skid to the trans- 
former and back to the negative bus. A 
quarter wavelength line with a natural 
frequency of 1 mc would be approxi- 
mately 250 feet long, which would appear 
to be more extensive than our circuit. 
However, because of the shape of the 
circuit—an irregular, tortuous loop— 
such a comparison is probably a poor 
approximation. The transformer itself 
oscillates at around 165 kc; its ampli- 
tude is about 50% of the high-frequency 
constituent. A yet lower frequency of 
about 12.5 ke is generated by the inter- 
phase transformer; this component also 
has about 50% amplitude, but because of 
its low frequency it persists much longer, 
and for this reason it may be potentially 
more dangerous. Very little of the 
voltage transient reaches other skids on 
the potline; the effect of the disturbance 
is felt only close to its source. 

Many other circuit conditions were 
investigated during the field tests; results 
followed very much the pattern of those 
described. For example, tests were made 
to simulate the effects of surge occurring 
at the time of an arc-back with both Y’s 
conducting. The only difference between 
this circuit arrangement and that just 
described was that a phase in the second 
Y was also conducting. 

The connection through the other Y 
was in parallel with all the other skids 
across the d-c bus. As this is already a 
comparatively low surge impedance, the 
extra parallel path could hardly be ex- 
pected to make any dramatic difference 
to the over-all picture; this indeed proved 
to be the case. It was interesting to 
observe slight differences in the responses 
according to which phase in the second 
Y was carrying current. These differ- 
ences arise because of the coupling, or 
lack of it (coupling is effected through the 
common primary winding), between the 


conducting phase and the arc-back 
phase. 

Two other test series are worthy of 
comment. The first showed the con- 
sequences of a surge after the primary a-c 
line had been opened, the second recorded 
events consequent upon a current surge 
entering a rectifier skid from its unit bus. 
The first of these conditions could arise 
when a potline is being de-energized; it 
would prevail if a surge occurred during 
the interval between tripping out the 
a-c circuit breakers and opening the 
cathode breakers and negative dis- 
connects. 


Under these circumstances we would. 
expect a difference in the low-frequency 
components of the response because the 
transformer presents a different surge 
impedance on open circuit than when 
short-circuited through the incoming 
supply.5® In general, the surge imped- 
ance is higher and the natural frequency 
lower for the open-circuit condition; 
this was reflected in the oscillograms. 
Fig. 9(B) is typical, showing that the 
natural frequency of the open-circuited 
transformer is approximately 50 ke in- 
stead of the 165 ke observed in Fig. 8, with 
the primary winding short-circuited. 
The high-frequency oscillations remain 
sensibly the same, so that the first peak is 
comparatively unaffected. In general, 
we may say that the effect of opening the 
a-c line is to cause the higher-voltage 
peaks to persist longer through the con- 
structive interference of the lower fre- 
quencies. 


For studying the behavior of a recti- 
fier skid when a surge enters from its unit 
bus, the skid concerned was isolated by 
opening the cathode breaker and the nega- 
tive disconnect. The stimulus was ap- 
plied and the response recorded at the 
unit bus, points 7 and 8. Two condi- 
tions were studied, normal operation and 
commutating. The transient excited 
when the circuit is in normal operation is 
shown in Fig. 9(C); the commutating 
condition was only slightly different be- 
cause of the two parallel paths through 
the transformer. The oscillograms in- 
dicate the generation of 67 volts/amp 
suppressed. 

Let us assume that there are three possi- 
ble ways in which a current surge can 
enter from the unit bus: 

1. It could be generated by a chop at 
the cathode breaker. 

2. It could come from a neighboring 
skid. 

3. It could be generated by the pot- 
line. 

The three possibilities will be considered 
in turn. 
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To generate an overvoltage of 1 kv the 
circuit breaker must chop 15 amp; 
to generate 10 kv, 150 amp, etc. Ample 
oscillograms of the behavior of the 
cathode breakers have been taken and 
there is absolutely no evidence that they 
chop in this manner. 

Regarding surges from neighboring 
skids, it may be said that all the skids in 
parallel would experience the surge; 
consequently the surge current would be 
only a small fraction of that in the surg- 
ing skid. This means that the over- 
voltage appearing in the surging skid 
would be much greater than that in all 
of the other skids. Thus, we need not 
look for trouble from this source. Finally, 
a surge of current from the potline would 
divide at the bus into parallel paths 
through all the connected skids in a more 
or less equitable manner. It would there- 
fore need a very considerable surge from 
the source to cause significant over- 
voltages at the unit bus and on the skids. 

In full operation only 1/14 of the surge 
ot suppression would be felt in any one 
skid or at any one transformer terminal. 
From Fig. 8(B) which shows the trans- 
former’s response to a chop, note that the 
peak voltage/amp suppressed is given 


by 


8 
Zoettective = 0.02 =40 volts/amp, approx 

To generate 35 kv, a figure that has 
actually been recorded, assuming the sys- 
tem is linear, would call for a suppression 
of 


35,000 


40 = 875 amp 


But this is the current per skid. If the 
suppression emanates from the potline, 
it must have a magnitude of 14X875= 
12,200 amp in the line itself. 

The very-high-frequency component 
of Fig. 7(B) does not contribute to the 
peak, being attenuated before that 
point is reached. However, the inter- 
mediate frequency is important, and to 
stimulate it significantly the current sup- 
pression should occur in a time less than 
half its period, ie., in 2 or 3 psec. It 
should be pointed out that because the 
low-frequency term makes the greatest 
contribution to the peak voltage at the 
transformer, a current suppression of the 
above magnitude taking as long as, say, 
30 ysec could produce a surge voltage of 
28 ky. 


Conclusions 


The authors believe that the method 
described is a very useful one for studying 
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the behavior of complex industrial power 
networks when they are disturbed by 
surges of current. Repetitive current 
surging already applied successfully to 
individual components’ is extended to 
an entire network. 

The following, more detailed, 
clusions apply specifically to the installa- 
tion tested. However, in a general way, 
they are applicable to other similar pot- 
line systems. The method, of course, 
could be equally well applied to many 
other kinds of industrial power networks. 

1. The overvoltages could be gen- 
erated on the system by abrupt changes 
in current in any branch in any skid. 
The effect would be felt only locally, 1.e., 
in the immediate Y to which the surging 
circuit was connected. Effects on the 
adjacent Y would be relatively slight, 
and on other skids, insignificant. 

2. The effective surge impedances are 
in the range of 30 to 65 ohms; as an 
upper limit these could give rise to 30 
to 65 volts/amp suppressed, so that to 
generate 35 kv a chop of at least 560 
amp would have to occur. 

As stated earlier, to determine the 
precise cause of surges was not the prime 
purpose of this investigation. Never- 
theless, deductions from the results 
strongly suggest the rectifiers as the 
probably source. We say this because 
of the absence of evidence that the anode 
or cathode breakers chop currents of the 
order of several hundred amp, and also 
because the possibility of dangerous 
surges originating in the potline is all but 
discounted by virtue of the very low surge 
impedance presented to the line by the 
many parallel skids and transformers. 


con- 


4. A surge in a rectifier would gen- 
erate several frequencies of voltage in 
local circuits. The very high frequency 
(1 mc) would only appear as a significant 
constituent when the change of current 
was extremely fast, 0.3 usec or faster. 
A somewhat lower voltage could appear 
across the transformer by current changes 
in 2 usec or less. When conduction exists 
in only one Y of a skid, the interphase 
transformer can present a compara- 
tively high surge impedance. This could 
give rise to voltage transients of the order 
of 35 volts/amp suppressed, and since 
the natural frequency of this component 
is low, the decay time of the current 
could be as long as 20 usec. Saturation 
of the interphase transformer will tend 
to limit this voltage. 


5. A surge in one rectifier tank with 
the adjacent tank (to which it is connected 
through the anode dividing reactor) con- 
ducting, would develop transient voltages 
in the local circuit comprising the two 
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tanks and the reactor only. Across the 
reactor surges of the order of 20 volts/ 
amp can appear. 

6. If the adjacent tank is not conduct-. 
ing when the tank surges, the highest volt-. 
age is generated across the tank. Pro-) 
ceeding toward the transformer, the 
amplitude of the surge diminishes. At 
the transformer terminals all the lower: 
frequencies are present together with! 
approximately 50% of the high fre- 
quency. 

7. It is usual practice nowadays to; 
protect the transformers by lightning; 
atresters from low-voltage terminals tox 
ground. Short-circuiting out the ar-: 
resters does not prevent the development 
of high voltages along the bus and at the; 
rectifier tank. The local loop formed by 
the bus work has sufficient surge imped-| 
ance to generate spikes of 50 volts/, 
amp. ; 

8. While, clearly, the use of arrester 
is wise, it would seem prudent to takes 
steps to reduce the disturbance at it 
point of origin, perhaps by placing aé 
capacitor across each rectifier tank. The 
surge impedance of the bus work is close ed 
to 100 ohms (estimated from the oscillo- 
grams making due allowance for theg 
apparent damping) and the natural fre+| 
quency is of the order of 1 mc. Theseg 
figures correspond to an equivalent L andd 
C of 1.6 phenrys and 0.0016 uf. A pro 
tective capacitor C of about 2.5 uf woul 
seem appropriate. 

This value, which may seem an 
suggested for the following reason: 1 
surge of 35 kv is reported to have beent 
measured before lightning arresters were 
installed. The points of measurement 
were the anode breaker and the negatives 
bus (points 9B and 15B) quite close ted 
where the arresters were subsequently 
located. If this was due to a current 
surge at a rectifier, we may expect that the 
surge voltage at the skid would be per+ 
haps 50% greater, say, 50 kv. Augmentt 
ing the stray capacitance by 2.5 uf will 
reduce the surge impedance by a factow 
of 40, so that the surge peak will be “ 
duced to 1.25 kv. 

If a higher level than this is acetate 
then, of course, less capacitance could be 
used. C’ should have in series with ii 
a resistance just sufficient to keep within 
reasonable bounds the discharge curren] 
from the capacitor when the rectifier firesy) 
Otherwise the connections to C’ should bel 
kept as short as possible. 
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HE COMMERCIAL introduction of 

the silicon-controlled rectifier in early 
1958 disclosed the possibilities of its use 
in many applications usually performed 
jby thyratrons, magnetic amplifiers, and 
jrotating amplifiers»? One such appli- 
cation, which so far has not been 
thoroughly explored, is the use of these 
ectifiers as the power amplifier in a d-c 
iservo system. This paper presents the 
esults of designing, building, and testing 
a servo system which incorporates the 
icontrolled rectifiers as the power ampli- 
(fier for the d-c servomotor. The final 
aim is to evaluate the performance of the 
ontrolled rectifier in this application. 

The silicon-controlled rectifier is basi- 
ally a p-n—p-n semiconductor with three 
ectifying junctions, in which regenerative 
switching is achieved through an ava- 
ilanche breakdown of the center p-n junc- 
tion. Breakdown occurs either by increas- 
ling the anode voltage to the breakover 
\point or by introducing a low-level gate to 


i . one 
Icathode current with the anode positive 
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but at less than the breakover voltage. 
Naturally, the latter method is used for 
power control applications. The rectifier 
regains its blocking state when the anode 
voltage drops down close to zero. This 
method of control is very similar to that 
used for thyratrons, except that the latter 
requires voltage control while the con- 
trolled rectifier uses a current control sig- 
nal. In addition, because of junction heat- 
ing effects and variation of firing current 
requirements with temperature, current 
pulses are more desirable than steady cur- 
rents as a means of firing the controlled 
rectifier. 

The controlled rectifier has the char- 
acteristics of small size, fast switching 
action, low voltage drop (about 1 volt), 
and high current rating, up to 16 
amp (amperes).*4 The thyratron, while 
possessing fast response and high current 
capacity, is bulky and has a relatively high 
voltage drop (10 volts or more). Thus, 
there is considerable incentive in investi- 
gating the performance of the controlled 
rectifier in applications currently per- 
formed by thyratrons, including the one 
investigated in this paper. Except for 
the firing control circuitry, much of the 
literature on thyratron control of d-c 
motors®® is applicable to the controlled 
rectifier, and offers a good starting point 
when the application of the controlled 
rectifier to a d-c servomotor is considered. 

The basic circuit selected for this in- 
vestigation is a half-wave bidirectional 
rectifier circuit. While producing less 
output current and more ripple than a full- 
wave circuit, the half-wave circuit is 
simpler, and adequate for d-c servo- 
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As mentioned previously, pulse firing 
is desirable for the silicon-controlled recti- 
fier. The basic control scheme used in 
the present investigation is pulse firing 
from unijunction transistors, which have 
the capability of performing as a relaxa- 
tion oscillator. This feature dovetails 
nicely into a simple pulse-firing scheme 
for the controlled rectifier. 

The circuitry and controls selected for 
the present system do not by any means 
represent an optimum, but rather an 
attempt to achieve a fairly simple system 
giving high servo performance. Thisisa 
very effective means of uncovering limita- 
tions, if any, imposed on servo perform- 
ance by the silicon-controlled rectifiers 
themselves. 


System Operation 


Fig. 1 is a schematic diagram of the 
over-all system. The servomotor is a 
small d-c torque motor rated at 115 volts, 
0.58 amp, 0.23 lb-ft (pound-foot). The 
motor is fed from the silicon-controlled 
rectifier circuit, which uses half-wave bi- 
directional control. The supply voltage 
is 115 volts, single-phase, 400 cycles. 
Two control signals determine the out- 
put of the rectifier circuit. One is the 
400-cycle a-c positional error signal 
from the control transformer. The 
other is a d-c stabilizing signal from the 
output of the tachometer amplifier and 
stability network. The latter signal pro- 
vides velocity feedback through a high- 
pass filter. The source of this signal is a 
d-c tachometer generator directly coupled 
to the torque motor. A control trans- 


IP,115 V 
400 ~ 22VDC DC TORQUE 
MOTOR 
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Fig. 2. Silicon-controlled rectifier circuit: power and control 


mitter, identical to the control trans- 
former, serves as the position reference 
for the servo system. 

The heart of the system is the silicon- 
controlled rectifier circuit, a circuit dia- 
gram of which is shown in Fig. 2. The 
rectifiers and motor are connected in a 
half-wave bidirectional circuit supplied 
at A and B by 115 volts, single-phase, 400 
cycles. Power diodes D1 and D2 are used 
to permit a common cathode connection 
of the silicon-controlled rectifiers. This 
permits some simplification of the control 
circuitry. A 22-volt d-c supply is fed 
across terminals Cand D. This provides 
power for the transistor circuitry. 

The firing of silicon-controlled recti- 
fiers SCR1 and SCR2 is initiated by posi- 
tive current pulses from unijunction tran- 
sistors Ul and U2 respectively. (For uni- 
junction transistor characteristics, see 
references 7 and 8.) The range of firing 
control for each rectifier is from 10 to 180 
degrees. The relative firing angles of 
SCR1 and SCR2 determine the magnitude 
and direction of the average motor current 
and hence torque. For zero torque, both 
SCR1 and SCR2 are prevented from firing 
until late (160 degrees) in their respective 
half-cycles. Some conduction of both 
rectifiers is desirable at standstill in order 
to obtain greater sensitivity at null. 

The control of the firing angle for each 
silicon-controlled rectifier is accomplished 
as follows: Consider the positive half- 
cycle of supply voltage (A positive with 
respect to B). Assume that a positive 
control voltage e, (made up of positional 
signal égc, stabilizing signal eg,, and a d-c 
bias signal) exists at the base of transistor 
Q1A. The collector current from tran- 
sistor Q1A, and hence from Q1B as well, 
will be approximately proportional to ¢. 
Thus the rate of build-up of the voltage 
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€c across condenser C1, during the positive 
half-cycle, is instantaneously proportional 
to the control voltage ¢;. When the volt- 
age across Cl builds up to 12 volts, Ul 
discharges C1 into the gate of SCRI, 
causing SCR1 to conduct. 

When SCRI1 fires, C1 is short-circuited 
by SCR1 via diode D1A. During the 
negative half-cycle, e-: remains clamped at 
about zero volts via leakage resistor R4 
and diodes Dl and D1A. When terminal 
A again becomes positive, the condenser 
voltage é can build up again in accord- 
ance with the net control voltage e;. The 
net result is that the firing angle of SCR1 
is inversely proportional to the average 
value of control voltage e, up to the point 
of firing, provided e is positive during the 
positive half-cycle. For e, negative, no 
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Fig. 3. Rectifier circuit waveforms 
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build-up would occur. With control 
signals eae and eae equal to zero, & has a 
d-c bias to provide a charging rate into 
C1 such as to fire SCR1 at about 160 
degrees. ; 
From the symmetry of the circuit it 
can be seen that whatever applies to con- 
trol voltage e¢, during the positive half- 
cycle also applies to control voltage e, dur- 
ing the negative half-cycle. Thus the 
firing angle of SCR2 is inversely propor- 
tional to the average of é2, up to the point 
of firing. Control voltage e, also has a 
d-c bias that is sufficient by itself to fire 
SCR2 at about 160 degrees. by 
Fig. 3 shows some circuit waveforms 
for various combinations of control signals 
€ac and ége. For the condition of zero 
control signals, as shown in Fig. 3(A), 
both SCR1 and SCR2 fire at the same 
angle (160 degrees) so that the net torque 
developed is zero. Fig. 3(B) illustrates 
the case of a forward (in-phase) error 
signal éac, with ea, and back electromotive 
force (emf) ep equal to zero. During t 
positive half-cycle, when condenser Cl 
is unclamped, eac is positive and causes 
condenser voltage e,:to build up at a faster 
rate than the zero-signal rate. This re- 
sults in an advancement of the firing me | 
for SCR1. During the negative half- 
cycle, when condenser C2 is unclamp 
€ac is negative and causes the firing of 
SCR2 to be eliminated. Thus a forwar 
torque is developed, tending to reduce the 
positional error signal ég. In Fig. = 


a reverse stabilizing signal eae is shown 
added. Since a negative eg tends to de- 
crease control voltage ¢, while increasing 
é, it has an effect opposite to that of the 
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ig. 4. Tachometer amplifier and stabilizing 
network 


orward error signal eg. For the case 
hown, éac has the larger effect and ad- 
rances the firing of SCR1, though not as 
much as before. Thus the forward 
orque is reduced. 

Fig. 3(D) illustrates the effect of motor 
pack emf. The same magnitude of con- 
ol signals eg, and ege as in Fig. 3(C) is 
ssumed. The presence of a negative 
back emf (forward motor rotation), delays 
he beginning of positive anode voltage 
n SCR1 while advancing it for SCR2. 
ince the firing angles are keyed with 
espect to the beginning of positive anode 
roltage, the firing of SCR1 is delayed while 
hat of SCR2 is advanced. In addition, 
he samme back emf cuts off SCR1 earlier 
vhile allowing SCR2 to conduct longer. 
astly, the instantaneous net voltage 
vailable to drive current through the 
motor armature resistance is less during 
HCR1 conduction and more during SCR2 
onduction. The result is that SCR1 con- 


flucts less average current, while SCR2 


The circuit diagram for the tachometer 
implifier and stabilizing network is shown 
Fig. 4. The bias potentiometer is 
djusted to obtair about 10 volts across 
ansistor output resistor R5 thus per- 


lbontrol circuit, forms a high-pass net- 
\work of the form 7,/(1+Ts) where T is 
approximately 0.03 sec (second). Con- 
Hlensers C4, C5, and C6 serve to minimize 
L00-cycle noise pickup as well as 400- 
bycle feedback from the torque motor 
ria the direct-coupled d-c tachometer 


| The locked rotor torque characteristics 
fan be readily derived on the basis of 
lhegligible armature inductance and recti- 
ier voltage drop, with the zero-signal 
tonduction angle ¢) as a parameter. 


The resultant characteristics show a 
similarity for both a-c and d-c control 
signals up to a torque output equal to 
one half of maximum. This is expected 
since the rectifier output depends on the 


. average of the control signal over the 


firing angle. The torque characteristic 
for ¢>=160 degrees is about the best as 
far as linearity and smoothness of con- 
trol are concerned; this is the angle 
selected as the zero signal conduction 
angle. For ¢)=160 degrees, the torque 
curve can be considered linear up to about 
40% of T maximum with saturation be- 
ginning thereafter. Locked rotor current 
measurements confirm this. 

The speed-torque characteristics were 
measured rather than derived, because of 
the added computing difficulty caused by 
the 250-ohm motor-shunting resistor. 
The resistor has the desirable effect of 
essentially reducing the circuit inductance. 
This in turn results in an increased locked 
rotor torque as well as a smoother transi- 
tion between positive and negative torque. 
However, the regenerative braking effect 
of the shunting resistor during the non- 
conducting part of the cycle modifies the 
speed-torque characteristics and causes a 
reduction in maximum speed. This is 
not objectionable for the present system, 
since the torque motor is not mechanically 
designed for high speeds. 

Measured speed-torque characteristics 
for low speeds and torque with ¢)=160 
degrees, are shown in Fig. 5. Tmax 
corresponds to maximum locked rotor 
torque and for this system is equal to 
0.07 pound-foot. The theoretical speed 
corresponding to a motor back emf equal 
to the peak of the supply voltage is de- 
fined as Q maximum. For this system 
Omax= 280. rad/sec (radians per sec). 
The a-c control signal a and d-c control 


signal 8 are normalized signals at the in- 
put to transistors Q1A and Q2A. They 
are defined as 


where ég is the peak of the a-c signal com- 
ponent, eg is the d-c signal component, 
and é is the d-c bias signal required to 
produce ¢e, all with respect to ground at 
the input to transistors Q1A and Q2A. 
This representation, together with the 
similarity of characteristics produced by 
small a-c and d-c control signals at low 
speed, permits the speed-torque curves to 
be shown as identical for both types of 
signals. The conversion of rectifier input 
signals eg, and ége to a and # respectively 
can be done from the measured relation- 
ships: 


€ac (rms)=0.07 @ volt 


€ac = 0.3 B volt 


At the low speeds and torques shown in 
Fig. 5, the characteristics are essentially 
straight, parallel lines. Also for —0.6< 
a (or 8) <0.6, the lines are approximately 
equally spaced. Thus, for this range of 
control signals and for —0.15<Q/Qmax 
<0.15, the system may be considered 
linear, 


Servo Analysis and Performance 


An approximate frequency response for 
the rectifier circuit (locked rotor) can be 
derived under the assumptions that the 
signal frequency is small compared to the 
supply frequency, that only one of the 


Tma x 


Fig. 5. Measured speed torque curves, low speeds; a-c control signal a, or d-c signal 8, as a 
parameter 


Cantor—Silicon-Controlled Rectifiers in a D-C Servo System 


SB etodees 
GR(CALC,) = € 1600 


|Ga| (DB) 


db 


°o 


[se (DEGREES) 


|Gg(CALC)|= |G,(MEAS)| = aN 


RECTIFIER 
CIRCUIT 


I I 
fo} 


Gp IS NORMALIZED WITH 
RESPECT TO ZERO FREQUENCY 
OUTPUT CURRENT Io 


| ee | 
aay 7 10 i00 1000 
Ww (RAD./SEC.) 
Fig. 6. Rectifier circuit frequency response 
RECTIFIER 
CON. TRANSF. AC GAIN RECTIFIER 
& GEARING CONTROL & MOTOR 


40 is 


hn 
ic} 
am) 


fe) 


|¢,| (oEcIBELS) —~ 
8 


67 
G, (CALCULATED) = —-—s 
S(l+g3) 
G, (MEASURED) 


-80 


(DEGREES) —>— 


10 
W(RAD/SEC.) 


Fig. 7. Uncompensated open-loop frequency response 


rectifiers is conducting during any one 
cycle of supply frequency, and that this 
conduction is small. Thus as an approxi- 
mation, we may consider the conduction 
to be a current pulse at the end of the half- 
cycle over which the input control signal 
is averaged. The final result is 


tL w 
Gua Fae aul T /4 


where T is the period of the supply fre- 
quency, J; is the fundamental component 
(signal frequency w-¢) of the rectifier circuit 
output current, and J) is the average out- 
put current at zero signal frequency (d-c). 
The same result is obtained for larger 
rectifier outputs that are still proportional 
to the input signals, if it is assumed that 
the conduction area can be approximately 
replaced by a current impulse of equal 
area located at the center of conduction. 

The fundamental J; at signal frequency 
w- is, strictly speaking, not a fundamental 
except when the supply frequency ws is an 
integral multiple of w-. However, even 
at in-between frequencies for w,, the 


10 


predominant harmonic will be the one at 
signal frequency we, as long as we is small 
compared tows. This major component is 
referred to as the fundamental [;, to 
which Gg=e’"*/* applies. The 1/4- 
cycle delay of this approximate frequency 
response is due to the averaging used in 
the control scheme and to the character- 
istics of the half-wave circuit. It should 
be noted that this delay does not apply 
to sudden inputs such as a step function, 
where the delay time can be as much as 
one cycle of the supply frequency. 

The calculated frequency response of 
Gr.is shown in Fig. 6. Also shown is the 
measured frequency response of the recti- 
fier circuit. The amplitude response is 
the same for both, namely unity. The 
phase response of Gz (calculated) and Gz 
(measured) shows good correlation, con- 
sidering the approximations used in de- 
riving Gp (calculated). In any case the 
phase lag of the entire rectifier circuit is 
very small up to 50 eps (cycles per sec) 
which is one eighth of the supply fre- 
quency. 
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Assuming that the motor load is pure 
inertia, that the armature inductance is 
negligible, and that the lag of the rectifier 
circuit is negligible, the uncompensated 
system open-loop response for small sig- 
nals and speeds can be obtained from the: 
speed-torque curves of Fig. 5. This 
proves to be 7 


Me 
~ a(or B) Ss . 
Los q 
( ‘a i 


The amplitude and phase response of 
this calculated G, are shown in Fig. 7} 
also shown is the measured frequency 
response. The differences in amplitude 
and phase are accounted for by the 
coulomb friction present in the system. 
The measured response is a further veri- 
fication of the negligible phase lag of the 
rectifier circuit itself, at frequencies up 
50 eps. 

An analytic approximation of the 
measured Gy, is 


40 


An apineedeiaa tion of the compensated 
open-loop response G is 


G =G for G,A,<1 


1 
G mE for G;Hi>1 


ee mnie i ee ep ante ation as 


Tn order that the compensated open-loop 
response may have a —20-db (decibels) / 
decade slope from about w=33 to #=500° 
rad/sec, the minor loop feedback A, is 
selected to be 
. 


0.0255? 


An asymptotic approximation of the 
compensated open-loop response G is 
shown in Fig. 8. It is desired that the 
complete open-loop response KG cross 
the zero db line at about w=350 rad/sec, 
This should result in a system of high 
bandwidth and good stability. There- 
fore = is adjusted to a value of 300, or 
Ko= = ‘ : 

The ign closed-loop frequency 
response of the system for small signals is 
shown in Fig.9. The response is flat to 43 
eps, reaching a peak of 1.4 (3 db) at 55 
cps, corresponding to the 90-degree lag 
point. The +3-db bandwidth is about 
65 cps. This response correlates fairly 
well with the previous open-loop meas- 
urements and calculations. 

Typical transient response measure- 
ments are shown in Fig. 10. The 1/4 
degree step function represents operation 


i 
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Fig. 9. Clcsed-loop frequency response, measured 


within the approximate linear range of the 
rectifier circuit. The 5-millisecond rise 
time is 0.33 of the period of the 65-cps 
bandwidth, which is reasonable. The 7- 
degree step function represents torque 
} saturation through most of the response. 
This signal does not produce velocity 
saturation. A step function of 200 de- 
grees or larger would be required to reach 
velocity saturation and this is not feasible 
with the single high-speed synchro in the 
existing system. However, other tests 
indicate that the maximum slewing 
velocity would be about 60 rad/sec. It 
should be noted that the transient re- 
sponse measurements do not include an 
initial delay time, which for the half-wave 
circuit can be as much as one cycle of the 
400-cycle supply frequency. 

'The measured velocity constant of the 
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system is 10,000 sec™4. This checks 


closely with the velocity constant of 12,- 
000 sec obtained from the open-loop re- 
sponse of Fig. 8 (with K =300). 

The torque constant of the system, with 
respect to the motor shaft, is 15-lb-ft/ 


rad. The coulomb friction in the system | 


is about 0.015 lb-ft. Thus the theoretical 
motor position error due to frictionis 1 X 
10% rad or about 0.06 degree. The ac- 
tual position error is much less than this 
due to the 400-cycle dither inherent in the 
system. The measured friction error is 
about 0.01 degree. 


Conclusions 


The use of silicon-controlled rectifiers as 
the supply for a d-c servomotor permits 
the design and realization of a servo sys- 
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Fig. 10. System transient response 


tem of very high performance. As far 
as servo performance is concerned, there 
is no time lag associated with it other 
than that inherent in the control scheme 
and half-wave rectifier circuit. It is 
possible, using these rectifiers, to achieve 
a servo bandwidth up to one sixth of the 
supply frequency, and in all probability 
a bandwidth up to one fourth of the 
supply frequency is feasible. 

The experimental system is designed 
around a small d-c torque motor. The 
same basic circtit, however, could be used 
for d-c servo motors up to about 1/2 
horsepower, assuming that the motor it- 
self had suitable servo characteristics. 
About the only changes required would be 
an increase in the heat sink capacity for 
the silicon-controlled rectifiers and an 
increase in the rating of the power diodes, 
such as to match the current requirements 
of the motor. The control] circuit could 
remain very much the same. 

The pulse-firing control scheme used 
gives very fast, reliable performance. 
During many hours of operation,there was 
no evidence of misfiring. The breakover 
voltage of the unijunction transistors re- 
mains stable at 12 volts when used with 
a constant 22-volt d-c supply. In the 
case of both the d-c and carrier-modulated 
control signals, the variable rate-charging 
circuit used gives a sharp intersection of 
the condenser voltage wave with the 
breakover voltage, thus permitting pre- 
cise firing. The flexibility afforded by 
the use of either d-c or carrier-modulated 
control signals is another advantage of 
the circuit. 

The main disadvantage of the firing 
control scheme is the measured 30-degree 
variation in zero signal firing angle during 
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circuit warm-up, due to variations in 
germanium transistor leakage with tem- 
perature. This produces a change in 
rectifier circuit gain, but very little drift 
since the balanced circuitry minimizes 
rectifier output drift. The variation in 
zero-signal firing angle can be minimized 
by the use of a lower transistor amplifier 
gain or by the use of silicon transistors. 
The d-c servo performance obtainable 
with silicon-controlled rectifiers is at least 
equal to that obtainable with thyratrons. 
This together with the much smaller size 
and lower voltage drop of the controlled 


rectifier, makes it superior to the thy- 
ratron for use as the power amplifier in 
a d-c servo system. 
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Variable 3-Phase Reversible Voltage 
Sources With Use of Saturable 


Reactors 


M. BOLT 


NONMEMBER AIEE 


Synopsis: This paper describes the circuit 
arrangements, with appropriate analyses, 
by which two or three saturable reactors 
may be used to supply a 3-phase load. 
The output voltages of the devices are 
constant-frequency 3-phase, almost sinus- 
oidal sources which may be smoothly 
and continuously varied from an almost 
balanced positive sequence to an almost 
balanced negative sequence by variation of 
the saturable-reactor control currents. 


EVERAL recent publications!—® de- 

scribe investigations in which the 
polyphase wound-rotor induction motor 
has been used as the basis of stepless 
a-c variable-speed drives. Each of these 
schemes is characterized by the use of 
some form of static magnetic power mod- 
ulator which supplies an induction motor. 
References 1 and 2 describe schemes in 
which saturable reactors are used in series 
with individual machine windings, while 
the scheme of reference 3 involves the use 
of static magnetic frequency doublers. 
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L. R. Foote, in reference 4, has dem- 
onstrated the novel idea of using 
saturable reactors as interphase magnetic 


‘couplings in order to produce certain 


desirable voltage unbalances for speed 
control purposes. 

This paper demonstrates some extended 
possibilities by which the interphase cou- 
pling of saturable reactors can produce 
variable, controllable, 3-phase voltages 
suitable for the control of polyphase in- 
duction motors in situations which re- 
quire a simultaneous application of posi- 
tive- and negative-sequence power. 


The Basic Saturable Reactor 
Connection 


Fig. 1 shows the basic saturable reactor 
interphase connection, classified for con- 
venience as connection 1. Two power 
windings, having the turns ratio p to 1 
are connected in series with phases B and 
C respectively, with the magnetic polarity 
as shown by the dots. 

To facilitate the analysis the reactor is 
assutned to be ideal, having no losses and 
no leakage flux. When this ideal mag- 
netic coupling has unity turns ratio, the 
self-inductances of each power winding and 
the mutual inductance between them are 
equal. The mutual reactance is sym- 
bolized by X in Fig. 1 and is a function of 
the reactor control current. 
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It is further assumed that the reactor 
core has an idealized magnetization char- 
acteristic of the type shown in Fig. 2. 
X is then variable from zero at saturation 
to infinity in the unsaturated region = 
has different finite values corresponding 
to every point on the curved transitional 
region of the magnetization character- 
istic. ! 

These assumptions permit the circuit 
of Fig. 1 to be analyzed by con- 
ventional steady-state sinusoidal. cit- 
cuit analysis techniques. The validity 
of the assumptions is discussed later. _ 

An outline of the analysis of connection 


Fig. 1. 


Basic saturable reactor 
connection 1 


5 
a 


Fig. 2. Idealized magnetization characteristic 
of saturable reactor cores | 
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Limiting load voltage vectors for 
basic connection 1 


1 is given in Appendix I for the case of a 
balanced load, Z,=R,;+jX, ohms per 
phase, and a unity turns ratio reactor. 

Equations 10 and 11 from Appendix I 
are reproduced here for convenience; let 


an Soe 10 
aos, Cy) 
=e: Vi Os 

Vee= ¥24| cos oe sin «le (11) 


The voltage vector diagrams for con- 
Mection 1 are shown in Figs. 3 and 4. 
With zero direct-control current the re- 
actor core has infinite permeability and 
tends to infinity. The corresponding 
load voltages are denoted by a, b(©), and 
ic(©) in Fig. 3. With maximum direct- 
icontrol current the reactor becomes 
katurated, having zero permeability and 

educing X to zero. Load voltages Van, 
Von, Ven, are then coincident with the 
applied voltages as denoted by the sym- 
bols a, b(0), c(0), in Fig. 3. These two 
extreme voltage conditions are independ- 
ent of the load impedance. 


. 


A,a 


Taking point C of the voltage vector 
triangle ABC as origin, and line CB as 
the x-axis, equation 11 may be expressed 
in terms of the Cartesian co-ordinates « 
and y of voltage Va. 


VoesX1,\? 
4R, 


_[ vea\* [ VowXr \7 
“( 4 ) +( ARE ) (12) 


Equation 12 describes a set of circles 
having their centers on the perpendicular 
through point V¢,/4 in triangle ABC. 
The perpendicular displacement y) of 
the circle center from line CB and the 
circle radius 7) are both functions of the 
load impedance. 

With a purely resistive load, the loci of 
phase voltages Von, Ve, as X varies from 
zero to infinity are semicircles on the 
diameters Vop/2 Vg¢/2 respectively, Fig. 
4(A) shows these loci with the phase 
voltage vectors at an arbitrary inter- 
mediate position of saturable reactor 
control current. Since, by equation 1, 
the load voltages must sum to zero and, 
by equation 6, Van is constant, vectors 
Von and Vé, must lie along adjacent sides 
of a parallelogram with diagonal NA. 

When the balanced load is partially 
reactive, the voltage loci Vj, and Vz,will 
follow the circular paths described by 
equation 12 and illustrated in Figs. 4(B) 
and (C). The limiting condition of an 
entirely reactive load produces loci along 
the line of vector voltage Voz. It is 
significant to note that the locations of 
the voltage loci described by equation 12 
are independent of the saturable reactor 
reactance. For a fixed load impedance, 
the value of thesaturable reactor reactance 
determines only the point of operation on 
a fixed locus. 


locus for an entirely 


locus for a 
inductive load 


resistive load N.n 


Let a different constraint now be 
applied in such a way that the magnitude 
only of the load impedance is kept con- 
stant, instead of both the magnitude and 
the phase angle. 


|Z,| =V/Ry?+X3? =constant (13) 


Combining equations 12 and 13 for a 
particular saturable reactor reactance X 
gives 


~ eye ta (Zee XY 
ae Sara Oma ye UA 
(15) 


The loci of constant-load amplitude 
described by equation 15 are a set of eccen- 
tric circles centered along the line of volt- 
age Voz. Since the radii and the loca- 
tions of the circle centers are functions of 
both the load impedance amplitude and of 
saturable reactor reactance X, it follows 
that, by controlling the value of X, a pre- 
scribed voltage unbalance can be sus- 
tained in the presence of load variation. 

The voltage vector loci of Figs. 3 and 4 
show that a load voltage phase inversion 
cannot be achieved by connection 1 for 
any R-L (resistance-inductance) load 
combination. The use of nonunity turns 
ratios with this connection results in 
different voltage loci but these show no 
advantage over the unity turns ratio per- 
formance described above. 

The principle and analytical technique 
used for connection 1 is now extended to 
more complicated saturable reactor con- 
nections. 


Reactor Connections Giving 
Reversible Output Voltages 


One of the axiomatic requirements of 
a variable-speed drive is that the driving 


A 


locus for a 
resistive load 


A 
(A) 
Fig. 4. Load voltage loci for basic connection 1 
A—Resistive load B—Resistance-inductance load C—Resistance-capacitance load 
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Fig. 5. 


A—Connection 2 


motor be reversible. Now, a polyphase 
induction motor is reversed when the 
phase sequence of its applied voltages is 
reversed. Therefore the output voltages 
of any magnetic power modulator which 
is used for driving an induction motor 
must, as a basic requirement, be capable 
of controlled phase inversion. 

Several interphase combinations of two 
ot more saturable reactors may be used 
to produce controllable, unbalanced volt- 
ages with the additional property of con- 
trolled phase inversion. For convenience 
these are classified here into two groups 
according to whether the interconnection 
involves two of the phase conductors or 
all three. 


REACTOR INTERCONNECTIONS BETWEEN 
Two PHASES 


Figs. 5(A), (B), and (C) show the 
schematic circuit arrangements of con- 
nections 2, 3, and 4 respectively. If 
ideal series-connected saturable reactors 
are used as the interconnection imped- 
ances Z; and Z, in Fig. 5(A) and (B), it is 
found that each of the connections 2, 3, 
and 4 has identical theoretical load voltage 


A,a 


Fig. 6. Load voltage loci for connection 2. 
with purely resistive load, when Z, is a series 
saturable reactor 
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Schematics of saturable-reactor interphase connections 


B—Connection 3 


loci for the same fixed load. The loca- 
tions of the load voltage vectors on these 
fixed loci are determined only by the 
various saturable reactor control cur- 
rents. 

The three connections 2, 3, and 4 have 
the additional common feature that the 
load phase voltage Van is at all times 
coincident with applied phase voltage 
Vuv. No. displacement of the load 
neutral point occurs, irrespective of the 
type of loading or of the nature of inter- 
connection impedances 2, and Zp. 

An outline analysis of connection 2, 
Fig. 5(A), is given in Appendix II for the 
general load impedance 7,=R,+jX, 
and the general interconnection imped- 
ance Z,=R,+jX1, used with a unity 
turns ratio reactor. Analyzed in terms 
of all three variables Z,, Z, and sattirable 
reactor reactance X, the load voltage loci 
are fourth power functions of forbidding 
complexity. When, however, a con- 
straint is placed in turn on each of the 
controlling variables Z, and .X, the load 
voltage loci are found to be circles similar 
to the circular loci of basic connection 1. 

The performance of connection 2, Fig. 
5(A), with reactance X infinite, is sum- 
marized by equations 22, 24, and 25, 
reproduced here for convenience. 


tan =| Xs— Ry ETOH / 
Rr+2Ri 
Xp +2Xy 
R,+X, ——— 
l Late ey (22) 
With R;=0 
Ver Xr 
Vigo SDs 
o= [cos R, sin al (24) 
With X,=0 
_Vce Ri 
Voo= 5 E oie a sin a lp (25) 


When the interconnection impedance 
Z, in Fig. 5(A) is an ideal saturable re- 
actor, Ri=0 and the voltage loci are 
given by equation 24. With a resistive 
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C—Symmetrical connection 4 


Ce i ee A 


load, X,=0 and the negative part of 
equation 24 disappear. Voltage Vc» then 
has a semicircular locus on the diameter 
Vop/2, the location of the vector being 
determined by saturable reactor react= 
ance X4. i 
With X, infinite, connection 2 becomes 
identical with the basic interphase con- 
nection 1. Points € and 6 in triangle ABC 
then traverse circular loci from their posi- 
tive sequence locations C and B respec 
tively when X=0 to coincidence 
Vop/2 and when X= ©, as shown in Fig. 
3. Let reactor X now be maintained at 
infinity while interconnection reactor 4 
goes into saturation. Load voltages 


B 
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Fig. 7. Limiting load voltage conditions 
for connection 2 i 

A=xX20 eae : 

B—X= ©, X,=0 i 


Fig. 8. Measured load voltage loci for 
symmetrical connection 4 with resistive load 


Von and Vz, then proceed from co- 
incidence at V,;/2 around semicircle loci 
to coincide with applied voltages Voy, 
V py respectively, thus forming a balanced 
negative-sequence set. The loci of these 
load voltages are shown in Fig. 6 with 
the two extreme conditions also shown in 
Big, 7. 

When this control sequence is applied 
to nonresistive loads, the voltage loci 
still perform excursions between the two 
limits of Fig. 7 but these loci, now 
described by the complete equation 24, are 
then of the type shown in Fig. 4. Con- 
nections 3 and 4 also produce voltage loci 
identical with those of Figs. 6 and 7, 
when appropriate control sequences are 
applied. 


MEASURED PERFORMANCE OF 
SYMMETRICAL CONNECTION 4 


‘Shown in Fig. 8 are the measured load 
voltage loci for a resistive load used with 
the symmetrical 2-reactor connection 4, 
from Fig. 5(C). 

Comparison of the measured and 
theoretical characteristics illustrates 
the degree of approximation involved 
in the idealizing of the saturable 
reactors. The unsaturated reactance of 
the two unity turns ratio reactors were 
large enough to give good agreement at 
the phase inversion point V¢,/2. It is 
seen that the measured characteristics 
are fairly symmetrical about chords 
Vac/2 and Voz/2. Equations 24 and 25 
therefore suggest that the saturable re- 
actors are, as might be expected, behaving 
like R-L circuits of approximately con- 
stant power factor rather than as lossless 
inductors. 

A minimum voltage unbalance is 
evident at each extreme control position 
in Fig. 8 owing to the nonzero saturated 
impedance of the reactors. This limit 
on the degree of voltage balance attain- 
lable may be minimized by using core 
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materials of minimum saturated perme- 
ability. It may be almost eliminated 
by using fixed capacitors in series 
with each power winding, designed to 
resonate at the maximum reactor control 
curvant setting. However, the possi- 
bility of a complete balance of applied 
voltages, though desirable, is not essen- 
tial for the induction motor in a variable 
speed drive. In many cases, because of 
the additional cost and bulk of capacitive 
compensation to eliminate a small per- 
centage unbalance, it may be considered 
an unnecessary refinement. 

The maximum voltage across the re- 
actor power windings is equal to the ap- 
plied line-to-line voltage for any R-L 
load. In Fig. 6, for instance, voltage 
Vee, Vav, Vaz and V¢ each has its maxi- 
mum values when coincident with vectors 
Vez or Vgg. A capacitive load would 
introduce the possibility of excessive 
steady-state reactor voltages as well as 
the undesirable possibility of subhar- 
monic ferroresonances. 

A measured harmonic analysis of load 
voltage waveforms showed that the de- 
gree of distortion was consistent with the 
characteristics of Fig. 8. With both re- 
actors unsaturated there was only 2% 
fifth harmonic, while at the corners B and 
C of triangle ABC, the most balanced 
condition, the waveforms contained 15% 
third harmonic and 5% fifth harmonic. 
The waveform oscillograms in this latter 
condition showed barely perceptible de- 
partures from the sinusoidal. At values 
of reactor control currents corresponding 
to the peaks of the semicircles in Fig. 8, 
a maximum harmonic: content of 40% 
third harmonic was measured, together 
with much smaller quantities of other 
odd harmonics. 


Nal 


VAN 
B 
(A) 
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Fig. 9. 
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SYMMETRICAL COMPONENTS OF 

CONNECTIONS 2, 3, AND 4 

The torque-speed characteristic of a 
polyphase induction motor supplied by 
unbalanced applied voltages may con- 
veniently be considered as the algebraic 
sum of two separate characteristics be- 
cause of the positive-sequence and nega- 
tive-sequence components of applied 
voltage acting independently. For ac- 
curate and sensitive control of both hoist- 
ing and lowering loads it is desirable that 
the unbalanced voltages be continuously 
adjustable to give any desired combina- 
tion of positive-sequence and negative- 
sequence torque. The thermal design of 
the motor and the duty cycle of the drive 
will determine the limitation of permissi- 
ble current unbalance. 

The relative magnitudes of the two 
sequence voltages given by symmetrical 
connection 4 with a resistive load are 
shown in Fig. 9(A). The theoretical 
characteristic is found to consist of two 
elliptical sections intersecting at the 
phase inversion point. Also plotted in 
Fig. 9(A) is the symmetrical component 
characteristic corresponding to the meas- 
ured locus in Fig. 8. 

When a purely inductive load is used 
with symmetrical connection 4, the 
theoretical symmetrical component locus 
is the linear characteristic shown in Fig. 
9(B). The theoretical locus for any R-L 
load will therefore lie within the shaded 
area in this figure. 

Since connections 2, 3, and 4 have the 
same theoretical load voltage loci, their 
theoretical symmetrical component loci 
are also identical with the characteristics. 
of Fig. 9. Connection 3, Fig. 5(B), has. 
the disadvantage of requiring three 
saturable reactors but gives highly sensi-~ 


2 
(B) 


Symmetrical component load voltages for connection 4 


A—Measured and theoretical loci with resistive load 
B—Theoretical boundary loci for resistance-inductance loads 
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Fig. 10. Saturable reactor interphase connections 


A—Cross connection 5 


tive control of the load voltages. Con- 
nections 2 and 4, shown in Figs. 5(A) and 
(C) respectively each require two reactors 
and appear to be equal in cost and bulk 
as well as identical in performance. L. 
R. Foote has described a scheme which 
is essentially that of connection 2 with a 
third reactor in series with line A. The 
symmetrical connection 4 is found, after 
considerable topological rearrangement, 
to be identical with a circuit suggested by 
Leonhard® and uses two identical reactors 
operating in identical modes. It has the 
advantage of resulting in a compact and 
elegant sinusoidal analysis. 


REACTOR INTERCONNECTIONS BETWEEN 
THREE PHASES 


Variable reversible 3-phase voltages 
may also be obtained by connecting two 
saturable reactors, as shown in Figs. 
10(A) and (B). These are classified as 
connections 5 and 6 and are referred to 
as the cross connection and the neutral 
connection respectively. 

An analysis of cross connection 5 for 
the general load Z,=R,+jX, with two 


Fig. 11. Load voltage loci for cross connection 5 with resistive load 


A—As X varies from zero to »; X;= 0 


16 


B—Neutral connection 6 


similar unity turns ratio reactors re- 
sults once more in complex fourth-power 
equations for the load voltage loci. To 
reduce the analysis to a manageable com- 
plexity, reactors X, and Xe, in Fig. 10(A) 
are considered ideal, and each in turn is 
constrained in value. The analysis, 
not included in the paper, gives the follow- 
ing results when X,;=©, X2=finite; 
then 


3Rr 


SUING ox are 26 
OED eT a 20) 
3V; xX 
Veo = ake (cos yee sin s) | 02 (27) 
R, — 
When Xo= © X,=finite, then 
3Ryz, 
we) el En ee ee 2 
DAES ES a uo) 
Xe 
Van = Vac{ jo sin a) | 0; (29) 
Ry, —_— 


It is evident that equations 11, 24, 27, 
and 29 are of identical form, differing 
only in the magnitude of the constant 
multiplier. In each case this multiplier 


defines the chord on which the circular 


B—Same as A for X2= 
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load voltage loci are based. For resistive 
loads the negative part of the equations 
disappears and the chord becomes a 
diameter. 

Fig. 11(A) shows the theoretical volt- 
age loci of cross connection 5 for a resistive 
load over the control cycle of reactor Xp 
while reactor X; is maintained at infinite 
impedance. The voltage vectors are 
shown at an arbitrary intermediate 
value of control current in reactor X9, 
With X,;= ~ X.=0, points band © coincide 
with points C and B respectively to give a 
balanced negative sequence. As X2 then 
increases, the load voltages 1 
circular loci, moving to points g when X= 
©, The load neutral point also describes 
a circular locus in which the line-to-line 
load voltage Vaz is at all times constant 
and equal to Vc. 

If reactance X2 is now maintained in- 
finite while reactor X; goes into satura- 


- tion, the voltages describe further circular 


paths from points g to points 7 in Fig 
11(B) which, for simplicity, shows onl 
the loci of points a and x. 


a 


A set of measured load voltage loci fo 
cross connection 5 are shown superposed 
on the complete theoretical voltage di 
gram in Fig. 11(C). Good agreement is 
evident only for the locus of phase voltage 
Vin. 

Voltages Van and Vo, exhibit a con 
siderable discrepancy from their respec. 
tive theoretical loci over the whole control 
cycle. It is evident that the saturable 
reactors are behaving nonideally and 
that the use of this connection as a power 
modulator in a close loop variable-spee¢ 
drive would require consideration in 
terms of measured rather than ideal 
theoretical performance. The use of 
cross connection 5 with nonresistive 
loads is governed by equations 27 and 29. 


x 


ee ee 


C—Measured and theoretical 
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Fig. 12. Load voltage symmetrical com- 
ponent loci for cross connection 5 with 
resistive load 


As with previous connections, the pres- 
ence of load reactance merely changes 
the radii and center displacements of the 
circular loci, which remain based on the 
chords shown in Fig. 11(C). 

It is found by experiment that the per- 
formance of neutral connection 6, shown 
in Fig. 10(B), is very similar to that of 
cross connection 5. 

Undesirably large neutral currents are 
necessary to achieve certain load voltage 
unbalances. This connection was not 
found to have any advantage over the 
performance of cross connection 5. 


SYMMETRICAL COMPONENT LOCI OF 
* Cross CONNECTION 5 


The measured and theoretical loci 
of the symmetrical component load volt- 
ages Vq, and V2, corresponding to Fig. 
11(C) are shown in Fig. 12. The lower 
theoretical locus, where X;= ©, is found 
to be elliptical, and identical to the 
corresponding portion of Fig. 9(A) for 
symmetrical connection 4. The upper 
theoretical portion of Fig. 12, where X2= 
© is seen to have a break point at 
rated applied voltage for the case of a 
resistive load R,. 

The considerable divergence between 
the measured and theoretical results in 
the voltage loci of Fig. 11(C) is reproduced 
in the symmetrical component loci of 
Fig 12. With a purely inductive load 
the theoretical voltage loci of Fig. 11(C) 
would become the chords of the present 
semicirclés. Fig. 13 shows the corre- 
sponding theoretical symmetrical com- 
ponent locus for cross connection 5 with 
an inductive load. The shaded section 
in Fig. 13 is thus the theoretical boundary 


sults consistent with the loci of Fig. 11(C). 
At points p for instance, where X,= ©, 
X2=0, load voltages Von, Ven each con- 
tained 17% third harmonic and 6% fifth 
harmonic while voltage Van contained 
only 1.0% of each. At points 7, on the 
other hand, voltages Van and Vz, con- 
tained 17% third and 5% fifth harmonic 
while voltage Vz, was very nearly a pure 
sinusoid. With both reactors unsat- 
urated, at points g, voltages Van and Von 
contained only 3% third harmonic. 

Comparison of the symmetrical com- 
ponent loci of Figs. 9 and 12, for connec- 
tions 4 and 5 respectively, indicate their 
relative usefulness as magnetic power 
modulators for driving polyphase induc- 
tion motors. The theoretical locus of 
cross connection 5, Fig. 12, shows 
promise, since positive-sequence voltages 
greater than the applied voltage are ob- 
tainable. 

Considering the two theoretical loci 
only, connection 5 appears to be more 
useful than the symmetrical connec- 
tion 4. It is evident, however, that the 
measured locus of symmetrical connec- 
tion 4, Fig. 8, approaches much more 
nearly its ideal theoretical performance 
than does the corresponding performance 
of cross connection 5, Fig. 11(B). Com- 
pared on a basis of measured perform- 
ance, the symmetrical connection 4 ap- 
pears to be the more useful. A final 
choice between the two connections would 
have to be made according to the quality 
of saturable reactors available and to the 
particular load voltages desired. 


Conclusions 


The paper has described five methods 
by which two saturable reactors may be 
interconnected between the lines of a 
balanced 3-phase supply to give control- 
lable, reversible 3-phase output voltages. 
Several of the connections appear to be 
directly applicable as power modulators 
for driving polyphase induction motors in 
stepless a-c variable-speed drives. 

The performance of the modulators 
can be predicted approximately by re- 
garding the saturable reactors as loss- 
less variable inductors and with the use 
of sinusoidal steady-state analysis. 


Appendix |. Outline Analysis 
of Basic Connection 1 


Fig. 13. Theoretical region of symmetrical 
components loci for R-L load with cross 
connection 5 


Combining equation 1 and equation 2 
yields 
Vap=1q(2Z,4jX)+1e(Z,+72X) (4) 


Combining equations 1 and 3: 
Vao= —To(Z,+j2X)—Ie2(Z,4+j2X) (5) 


Combining 4 and 5: 


2 
Vann =aTae Nee =constant (6) 
Combining 4, 5, and 6: 
Ven =1tZ1 
A Ven(2Z,+jX)+ Vea(Z1 +j2X) (7) 
3(Z,+j2x) 
Put 
Z,=Rr+jX1 
Vea = Ven(—1/2+j-V3/2) 
Vot = Ven— Ven (8) 
Combining 7 and 8: . 
ee 
c= a 
ON WS Ri2+(2X +X1) 
Ry 
tan“! ——~__ (9 
SS am orc 
Let 
Ry 
=———_ 10 
tan a aX +X, ( ) 


Combining 9 and 10: 


Vi = 2 cos pee sin « |le 
Hao R, os 


(11) 


Cartesian co-ordinates x and y of Vct 
are: 


L 
cso re ic analysis of th ares eat - e Yee [| Peas (12) 
A measured harmonic analysis of the Van <LoZp— Tol Zu-+4X) +10) X (2) ; ie 
load voltage waveforms for cross con- 
nection 5 with a resistive load gave re- Vac=1p(Z14j2X)—Ie(Z,+j2X) (3) Now, 
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|Z,| =VRi2+X2? (13) 


Combining 12 and 13: 
Vee |? Ver xs) 
| == A | “lr cE 


-[ MeelZsh 
| eRe 


Combining 13 and 14: 


2X* Vos |, 5 Veo ¥Zst1 
*—4x*—Z;3| > [4x*—|Z;|? 


15) 


Appendix Il. | Outline Analysis 
of Connection 2, Fig. 5 


(16) 
(17) 
(8) 
(19) 


Ig+Ip+le=0 

Vaz =loZ,—1y(Z1+jX) +1cejX 

Vac=ly( Zi +j2X)—IcZ,—Icj2X 

Vac =i X +10Z:—Ic(Zi1 +X) 
Combining equations 16, 17, and 18 yields 


2Vast Vac 
3 
=independently of 21 


Van=laZL= 
(20) 


If the value of X is now kept at infinity 


(21) 


(22) 


V/V R2+-X1? |p 


VJ (Rr +2R1)?+(X1+2X1)? 
(23) 


Voo= | Ves | 


(24) 


R 
Voy= V2 cos os sin a lle 
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Plant Identification in the Presence of 


Disturbances and Application to 
Digital Adaptive Systems 


P. JOSEPH 


STUDENT MEMBER AIEE 


HE NEED for an adaptive control 
system occurs when the parameters in- 
volved in the design of the system change 
in time. For small variations, conven- 
tional feedback configurations can some- 
times reduce the sensitivity of the output 
to these changes to acceptable values, 
but large variations of the parameters re- 
quire a more flexible system. The newer 
adaptive point of view in designing a 
control system can, in general, allow for 
changes in input signal characteristics, 
plant characteristics, and disturbance 
characteristics with a minimum amount of 
a priort knowledge. It appears that a 
general distinction has been made be- 
tween adaptive controls? and optimaliz- 
ing (optimizing) controls** although the 
characteristics are not always distinct. 
A recent note’ has suggested a possible 
classification according to whether there 
is an original feedback loop or not. 
Most methods of designing control 
system compensators assume that the 
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J. LEWIS 


ASSOCIATE MEMBER AIEE 


J. TOU 


MEMBER AIEE 


open-loop plant transfer function is 
known. If the open-loop plant transfer 
function could be determined, one of 
many compensator design methods could 
be mechanized to yield a plant-adaptive 
control system. A major problem in this 
approach to adaptive control, therefore, 
is the determination of the transfer func- 
tion: the identification problem. 

Any realistic treatment of the identifica- 
tion problem must include the effects of 
disturbances upon the accuracy of the 
identification. The presence of disturb- 
ances prevents an exact identification in a 
finite time and, therefore, a compromise 
between accuracy and speed of identifica- 
tion is necessary. 

In order adequately to describe a plant 
by a transfer function, the time variation 
of the plant parameters must be slow in 
relation to the significant system time 
constants. The method of identification 
used here actually assumes that the pa- 
rameters are constants between periods 


and the number of zeros of the transfer 


of adaptation. The relationship between 
rate of change of plant parameters and 
period of adaptation depends on many 
factors such as the identification accuracy 
required and the ratio of signal to output 
disturbance. For example, assume that 
a preselected identification accuracy of 
1% is chosen and that the signal/dis: 
turbance ratio is such that the time re- 
quired for this identification is 10 sampling 
periods. Ifa plant parameter is driftin g 
at a rate of less than 0.1% per sampling 
period, it can be assumed to be constant 
to within the identification accuracy 
during the identification interval. 

In carrying out the identification, it 1s 
assumed that both the number of poles 


function are known. If only an upper 
bound on each of these numbers is known, 
it is sufficient, although it is desirable to 
keep these bounds small since the com- 
putations required increase rapidly. 


Identification 


Consider the system of Fig. 1. G(g) is 
the open-loop plant transfer function and d 
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(2)= (1) 


V(z) is a disturbance and the output, 
(2)=V(z)+U(z). The input sequence 
R(z) and it is assumed that the input 
id the disturbance are independent. 
is assumed that the past values of m 
e known exactly since they are generated 
y the discrete filter which will be realized 
ly adigital computer. It is also assumed 
nat the past values of r and c¢ can be 
easured. Measurement errors in r will 
ot adversely affect the adaptation 
rocess and measurement errors in ¢ will 
ave an effect similar to that of w. 
herefore, they will not be considered, 
is assumed to be unmeasurable. 

If 


M(z) (2) 


mor, ares (3) 


where £ and 7 are known constants 
pr at least an upper bound of & and 7 are 
known. By using equations 2 and 3, 
2 set of simultaneous linear equations in 
he +7+1 unknowns, a, and };, can 
e derived. Any measurement error in ¢ 
an be reduced by the use of least squares 
smoothing. The success of the smooth- 
jng depends upon m being independent of 
the measurement error.’ 

) When U(z)+40, one might attempt to 
»stimate G(z) by G(z) where 


| _ Caz) 
| =a) (4) 
is pCa Os). U(z) 

(2) Uz)’ Mie) OF ais) (5) 


The use of least squares smoothing does 
not help in this case, since m is not in- 
dependent of ~. This is a result of the 
osed-loop nature of the system. 

As a starting point, an alternative to 
equation 1 is 


_ B27(2) 


©) Drm (2) 


(6) 


here ®,,,(z) is the correlation sequence 
between x and m defined by the relations 


i" 1 
brm(k) = fee ONL. > XM i+% (7) 


' G(2)= 


Fig. 1. System diagram 
R(z) 
and 
Pym(z) = ei dbzm(k)z * (8) 
k= —o 


and where x is any signal which is corre- 
lated with m. 

When U(z)#0, one might attempt to 
estimate G(z) by G(z) where 


A ®r0(z) 

Gere es (9) 
A _ Bro 2) + BrulZ) _ ®rx(Z) 
See A ERG. 


If x is chosen to be a signal that is not 
correlated with u, G(z)=G(z). The ob- 
vious choice for xis7.. Then 


(GB) es 
tO 


The solution of the identification prob- 
lem cannot, however, be left here. The 
evaluation of ®,,(z) and ®,,,(z) requires 
taking an infinite limit. This is not only 
impractical actually to implement, but 
would utilize information from the dis- 
tant past that would be out of date if 
the plant were time varying. If only 
recent information is to be utilized, an 
exact determination of G(z) is impossible. 
However, a good approximation is avail- 


(11) 


able. The following definitions are, 
therefore, made: 

I 
Orm(k, N,I)= >) remasx (12) 


t=I-—N 
foo} 


@mleN.I)= >) Om(k, N,D)s* (13) 
k=—0 
Reference 8 contains a proof of equa- 
tion 6, and a completely analogous argu- 
ment shows that 


G(z)= eee (14) 
G(z) will be estimated by G(z) where 
ate aT (15) 
If 

U(z)=0 

then 

G(z) =G(2) 

If 

U(z)40 


COMPENSATOR 
E(z) 


PLANT U(z) 


then 
sortase 
=G(s) ee as) 
Since 
aa Ra @ru(z, N, 1) =Sru(2) =0 


in the limit as NV goes to infinity 
G(z)=G(z) 


Before discussing the closeness of the 
estimation when JN is finite, equation 15 
will be manipulated into a form which can 
be used to solve for the 4,’s and 6,’s which 
determine G(z). G(z), the estimate of 
G(z), will be expressed in a form similar 
to equation 3 


£ . 
ajz 7 
¢@)—=—_ (17) 
14D) be? 
J=1 
& foo} 
Soe? SS adh Ne 
A j=0 k=—o 
G(z)= — = foo} 
1+) be DL Orm(b, N, Te 
j=1 k=—o 


Cross multiplying, letting k= —j—p, and 
equating coefficients of like powers of 2, 
yields — 


E 
= OO Ia Ps N, Dy 
j=0 


n 
>> dire —J—P, N, 1) =Orel—P, NT) (19) 


FS 


By selecting &+7+1 distinct values of 
p, &+n+1 linear simultaneous equations 
can be obtained from equation 19. These 
equations can then be solved for the 4,’s 
and the b,’s which determine G(z). 

Returning to the problem of the ac- 
curacy of the estimation of G(z) by G(z), 
or more precisely, the estimation of the 
a;’s and b;’s by the 4,’s and bys, it is found 
by examination of equations 14, 15, 16, 
and 19 that the a;’s and b;’s satisfy 


E 
4j9rm( —j—P, N, I)— 
j=0 


n 
> bir -J—, N, DT) =6re —B, ND) 
sa (20) 
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R(z) E(z) 


E 
SS 059 rm( —j—P, N, iN)es 


j=0 


7 
D> bylOrel Jb, N, 1)— 
j=1 : 

Oru —j—b, N, 1] 
=6,-(—), N, I) —6ru( —2, N, T) (21) 


Comparing equations 21 and 19 shows 
that the coefficients of by and by differ 
by Bru —j— p,N,1). Clearly 


El6ru(k, N, I)] =0 (22) 


and in Appendix I it is shown that 


=E{ [6ru(k, N, 1)]?} 
I iE 
=>) Dy brli-buali—f) (23) 
Ng 
and if u consists of white noise with 
variance o,”. 


oe = (N+ 1)$rr(0) on? 


Note that the square of the mean value of 
6re(k, N,1) is proportional to (V+1)? and, 
therefore, the percentage mean square 
difference between the coefficients of 
equations 19 and 21 decreases as 1/N-+1. 

Using a white noise disturbance, as an 
interesting special case, the following rela- 
tions are also developed in Appendix I. 
This development contains several ap- 
proximations. 

if 


(24) 


A?>K%6,,(0, N, I)ouw>) Ais? (25) 
aie 

then the lo error in the estimation of the 

a;’s (or bj’s) is 


Loerror = { E[(as—y)?] rcs (aj2-+1)7 
(26) 


where A is the determinant of the coeffi- 
cients of the simultaneous equations and 
where the A,,’s are the cofactors of the 
coefficients of the b;’s and K is a design 
parameter, _ 

The method of calculating the 0,m(k, 
N,J) in real time remains to be described. 
First, if the 6’s are to be evaluated at the 


DISCRETE 


2() 


ig. 2. tem diagram 
iE Fig. 2. Syste g 


C(z) 


nth sampling instant, examination of 
equation 12 shows that J<nandI<n—k 
since the values of the samples which will 
occur in the future are unknown. In 
order to utilize the most recent available 
information, the equality will be chosen 


IT=n—Rmax (27) 


where —kmax equals the most negative 
value of p that is selected. Suppose that 
the lower limit of the summation in equa- 


tion 12 is m) —Rmax then I—N=—FRmax 
or N=I—m+Fmax: 
N=n—nN (28) 


Note also that 


Orm(R, n—No0, n—Rmax) 
=Orm(R, n—1—m, n—1 —k)max+ 


1n—kmax™n+k—kmax (29) 


Starting with 
Orm(R, 0, No— 


Necro) = np —kmax™ ny) — kmax+k 


(30) 


the @’s can be calculated using equation 
29 which requires one addition and one 
multiplication, per 6, per sampling instant. 
As time, ”, goes on, JV increases. 

A result similar to equation 19 also 
can be obtained for continuous-data sys- 
tems and is shown in Appendix II. 


Application:to an Adaptive System 


The method of identifying the plant 
transfer function can be applied to adap- 
tive control in a number of ways. The 
control or command signal can be ob- 
tained as the output of a discrete com- 
pensator which is realized on a digital 
computer. 

The compensator must be changed in 
accordance with changes in the plant 
(or signal or disturbance) parameters 
and at least two general approaches to the 
adjustment of the compensator are possi- 
ble. In one approach, the compensator 
is redesigned each time the system 
“adapts.” The redesign is carried out for 
some preselected optimum type of com- 
pensation so that the performance 


Fig. 3 (left). 
System of the 
example 


Fig 4 (right). 


Discrete system 
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criterion is inherent in the compensator 
design. One alternative procedure is ta 
choose a performance criterion and cal- 
culate it from an identification of the 
closed-loop system. The compensator 
is then adjusted to give an extreme value 
of the performance criterion. This re- 
quires the mechanization of an op- 
timalizing loop which is inherently slower 
than the first method because it utilizes 
a searching technique. Either the re- 
design or the adjustment can be done 
periodically or according to some plan. 


The accuracy of the identification is 
determined from equations 25 and 26. 
A knowledge of the maximum tolerable I¢ 
errors in the a,’s and b,’s determine K by 
equation 26 and the plan is then to check 
equation 25 each sampling instant. When 
equation 25 is satisfied, the equations 
derived from equation 19 will be used to 
determine G(z). Since the expected value 
of the left-hand side of 25 is pro- 
portional to N+), while that of 
the right-hand side is proportional to 
N21, equation 25 is bound to be 
satisfied eventually unless lim A= 0. 


N-oO 


‘To avoid this, the values of used in 


equation 19 must be carefully selected to 
obtain a nonsingular set of equations. 
This is demonstrated in the example. 
After G(z) has been calculated and the 
compensator redesigned, the 6’s are cal- 
culated anew using equations 29 and 30 
with a ae m and the re 
repeats. 

If the maximum tolerable lo errors 
are not known a priort, the parameter 
K could be set by an optimalizing loop, 
or in the words of a recent paper, by i 
supervisor’s supervisor.’” | 


A large reduction in the number o! 
simultaneous equations that must be 
solved is possible if G(z) can be writter 
as G;(z) G,(z) where G;(z) is fixed ane 
well-known and G,(z) is varying and un: 
known. The system can then be re 
drawn as in Fig. 2 where Q(z)= CHO) 
M(z) and is easily calculated. 

Writing 


é 
Ds AjZ 
C\2)\= = 
1+) By? 


and similar to equation 19 
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Fig. 5. Time responses for system with changing plant-gain constant 


Initial Ag=10 
Rate of change of Ap=0.1/sample 
+ Adaptation instants 


& 
: : A j6rq(—j—P, iN f)— 
=0 


n’ 
>) Bite —ji—P, N, 1) =e —p, N, 1) 
j=1 
(32) 


G,(z) and hence G(z) can be determined 
by solving only &’++n’+1 equations, where 
&’<£ and n’<n. 

Although equations 32 and 19 utilize 
correlation functions with 7, it would 
also be satisfactory to correlate with a 
component of ry. It is usually preferable 
to use the nonconstant components of r 
only. 

If the disturbance does not occur at the 
output but elsewhere in the plant, the 
block diagram could still be manipulated 
linto the form of Fig. 2 where U(z) would 
now be an equivalent disturbance. In 
this instance U(z) would probably not be 
white noise. In many situations, such 
as position control with torque disturb- 
ance due to a steady wind or flow control 
with a disturbance due to leakage, the 
disturbance resembles a deterministic 
signal more than it does white noise. If 
this is the case, og? will be much larger 
than predicted by equation 24. How- 
ever, if a white noise test signal p is 
added to r and all correlations are per- 
formed with p instead of r, then 


a9? =(N +1) bpp(0)ou* 


This is because equation 23 is symmetrical 


(33) 
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© Input 


with respect to p and uw. Then all equa- 
tions and developments hold if 7 is re- 
placed by p. The introduction of p up- 
sets the system; however, it can be made 
very small if long periods between adapta- 
tions are tolerable. 


Example 


Consider as an example the system of 
Fig. 3. The completely discrete version 
of this system is shown in Fig. 4. 


Ayp+Aiz7 
Gos) Bie 
f'=y/=1 (34) 


There are three unknowns and an entire 
set of equations is available from equation 
32. It is necessary to select three of 
them. Consider the equation for p=0. 


Or(0, N, DT) Ao+6r¢( = 1,N,I)a1— 
Ore( —1, N, I)B, =6,-(0, N, T) (35) 


Suppose 7 consists of white noise, then 
due to the unit delay in the system 
Grq(0), Yra(—1), Pre (0), and Pre( — 1) equal 
zero. Thus if equation 35 is selected as an 
equation Lim A=0. Thus the equa- 


N-wo 
tions for =0 and similarly all p>0 are 
rejected. Since —max is the most nega- 
tive value of p selected and 


IT=n—kmax (27) 


it is desirable to choose the least nega- 
tive values of p possible in order to utilize 


X Output with adaptation 
A, Output without adaptation 


the most recent information. Thus p= 
—1, —2, —3 are selected. Then 


6r(1, ? N,I)Ay+6,,(0, N, I)A\— 
Ore(0, N, I)Bi=6rc(1, N, I) (36) 


Orq(2, N, T)Ay+6r(1, N, DAi— 
Ore(1, N, I)Bi=6;<(2, N,I) (37) 


679(3, N, 1)Ao+6r¢(2, N, I)Ai:— 
Ore( 2, N, 1B, =6r(3, N, T) (38) 


These equations are solved for Ao, 41, and 
Bi when equation 25 is satisfied. 

To illustrate the design of the compen- 
sator, assume deadbeat response to a step 
input is desired. If all the poles and 
zeros of G,(z) are inside the unit circle, 
i.e., |41/Ao|<1, |B,|<1, a straightforward 
application of the design procedure for 
deadbeat systems (reference 8, pp. 502- 
12) yields the following result: 


aotayz7} 


Ds) (39) 
where 

ay=1/A) 

a1=B,/Ao 

Bi=A1/Ao (40) 


Values for ao, a1, and 8; can be found sim- 
ilarly for the three other possible pole 
and zero configurations. After Ao, Ai, 
B, are determined from equations 36, 37, 
and 38, the compensator, D(z), is designed 
automatically in the computer with the 
aid of equation 40 and similar relations. 
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° ? 
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5 Percent 


02 
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K 


Fig. 6. Accuracy of gain (Ao) 


Computer Simulation 


The system of the example was simu- 
lated on a digital computer. An input of 
the form which is shown in Fig. 5 and a 
white Gaussian disturbance were used. 
Figs. 6, 7, and 8 show the experimental lo 
errors obtained in the determination of the 
gain and the zero and pole locations for 
various fixed values of K and signal— 
disturbance ratio. The values of lo 
errors predicted from equation 26 are 
shown by the straight lines. About ten 
data points were averaged to determine 
each point on the graphs and the trend 
of the errors to be inversely proportional 
to K and to be independent of signal— 
disturbance ratio is demonstrated. The 
ability of equation 26 to estimate the 
errors to within a factor of 2, justifies the 
approximations that were made. 


The period between adaptations, J, 
depends in a complicated manner on 
the shape of the input, on D(z), and on 
G(z). However, in order that the reader 
may form an estimate of the adaptation 
period, it may be noted that under the 
particular conditions used for the com- 
puter runs with K=100 and with a 
signal—disturbance ratio of 26 db (decibels) 
the adaptation time varied from 25 to 28 
sampling periods. If one notes that the 
left-hand side of equation 25 is propor- 
tional to (0,2N)***"*) and the right- 
hand side to o,?K2(o,2N)24*41, we 
see that the period between adapta- 
tions is proportional to K*o,2/a,2._ Since 
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500 - 


the errors are proportional to 1/K, the 
compromise between speed and accuracy 
of the identification is clear. 

An example of the performance of the 
adaptive system when the plant has a 
time varying parameter was obtained by 
allowing Ao in the transfer function 


Ao +A,zg71 


ied hose 


to increase at a steady rate of 0.1 per 


Fig. 7. 


Accuracy of zero (—Ai/Ao) : 
- 


sample. The initial value of Ap was 10 
which means that the parameter w. 


rectangular wave used previously w: 
again used as an input. The response 


the response of the system when t 
compensator was left at its original cor- 


is also shown. It can be seen that th 
adaptive system continued to have | 


ACCURACY OF POLE ¢B,) 


© 20 db Sig)... 


Error 


4 26 db Sign... 
a 40 db Sig, 


Percent 


Fig. 8. Accuracy 
of pole (—B;) 
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tisfactory performance whereas the 
pnadaptive system became  unsatis- 
ctory after about 20 sampling periods. 
he slight deviations from correct dead- 
eat response in the adaptive system are 
ie to the disturbance in the output 
hich causes errors in the identification 
5 well as actual disturbances to the sys- 
;no attempt was made to compensate 
> minimize disturbance effects in this 
ase. Some deviation from deadbeat 
sponse is also due to the fact that 
nce the parameters of the plant are con- 
antly changing and the compensator is 
Ways designed on the basis of past in- 
prmation, a small error in the compensa- 
br results. 


onclusions 


A method for identifying a linear plant 
1 the presence of an output disturbance 
as been presented. An approximate 
malysis has shown that the per-cent error 
identification can be made proportional 
D a parameter 1/K while the measure- 
nent interval is proportional to K?. 
imulation of a system utilizing this 
method of identification provided ex- 
erimental confirmation of these relation- 


Appendix |. Estimate of Error 
in Identification 


i I 
[6-u(k, N, nie] = rates | 


t=—-N 


=E[r]E[u)(N+1) (41) 


mee ry and uw are assumed independent: 
if either 7 or u contains no constant com- 
fonent, equation 41 equals zero. This 
issumption will be made. The variance 
f O;4(k, N, I) is 


f= 2( 10,.(k, N, I)}?) 


4 4 
= ( > ya restores) (42) 


Due to the independence of 7 and u 

I I 

f= >, >> enli-J)buli—f) (43) 
#=I-—N J=1-N 


f u consists of white noise, i.e., 


or 8X0 then equation 43 becomes 
@? =(N +1)¢r7(0) ou? 

ince 

[6,,(0, N, I)] =(N +1) ¢rr(0) 


1 estimate of ap? is 


(44) 
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oo’ =0;,(0, Ni I) ox? (45) 
Consider the set of equations 
(ori Qaione Bal pal 
=" (46) 
An... Ann Xn Yn 
By Cramer’s rule 
Yi O12. . . Min 
m=" nr. -Ann| _ D1 (47) 
QT eka vars ain A 
OTS 6 ie is lars Ann 


If the a;j’s in equation 46 are replaced 
by aij+(da)ij where (6a)Kay;, the errors 
introduced into x; are of interest here. 
Consider the error introduced into the 
denominator A. ; 


OA 


Oaij 


= Aj;=the cofactor of aij (48) 


The error introduced into A by (éaij;), is, 
therefore, A;;(6a)i;. If all of the (da);j’s 
are independent statistical quantities with 
standard deviation, o;;, the variance of the 
error introduced into A will be 


on? =E[(5A)?] =>) ye Aij?ois? 
ik id 


The assumption of the independence of the 
(6a);7’8 is not strictly applicable to the 
estimation of the errors in the solution of 
equation 19 but an approximation is all 
that is desired. 

Since 


dD; 
A+6A 


(49) 


is approximately equal to 
ID 1D} 


ran 


6A 
the error in x; due to 6A is 
6A 
be Ne 
Continuing the approximation assume 
oaA—op1 


then the approximate lo error in 4; is 


(Bl8m)}/2== ea(mt+l)/=0n (80) 
If 
A?<K*op? 
then 
oa<ae (mP+1)Y? (51) 
K 
where the condition can be rewritten 
(52) 


A2<K? ee ys Aig?ois? 


Appendix Il. Identification for 
Continuous Data Systems 


In order to measure the continuous 
transfer function, G,(s), of a plant, it will 
be assumed that G,(s) can be written as 


E 
ys ajs! 


Gaks) sods Ok Adee 


n 
1+) bys! 
f=) 


where £ and y are known. In writing 
equation 53, it is assumed that G,(s) con- 
tains no integrations. Equation 53 can 
be rewritten as 


(53) 


E 2 
ajs’_” 
G(s) == (54) 
s 14> bys 
j=1 
Since 
®,(s) 
Gis) == (55) 
(5) FAG) 
& 
ajs? "@,o(s)=s 1Br¢e(s)+ 
1=0 


n 
D5 bys! 7B re(s) (56) 


bei || 
Note that the inverse Fourier transform of 
s "®r¢(s ) 
is 
bran 7) 


where 


gx(t) 
t 


t 
= i} ane f q(t)dt* with k integrations 


foe} 


Thus the inverse transform of equation 56 
is 


é 


De agora —§(7) — 


=0 


n 
SD byhre q—9f T) = dren( py. (@sx/) 


by 3 | 

By selecting £+7-+1 distinct values of 7 
t+n-+1 equations can be obtained. This 
method requires 27 integrators to calculate 
the gx’s and cx’s plus (£-+-7-+1)? multipliers 
and integrators to perform the correlations. 
In addition, facilities to solve the simul- 
taneous equations must be provided. 
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Studies of Electrodynamic Forces 


Occurring at Electrical Contacts 


A. C. SNOWDON 


MEMBER AIEE 


NOWLEDGE of the magnitude of 
the electrodynamic force occurring at 
closed electric contacts carrying high cur- 
rent is important to the design of reliable 
contacting devices. Past work in this area 
by Holm! has shown that current passing 
through a contact surface is constricted 
to flow through very small contacting 
spots, and that this constricted current 
flow causes a repulsion force to act on the 
contacts. Further contributions to the 
theory and the calculation of this force 
have been made through studies at the 
Cutler-Hammer Laboratories in Mil- 
waukee, Wis.,2 the progress of which 
is reported in this paper. 

The following analysis derives equations 
to calculate the magnitude of the re- 
pulsion force just mentioned. It will be 
shown that the force depends on contact 
pressure, diameter, material hardness, and 
magnitude of current. The force is found 
to be independent of contact surface 
curvature, when contacts that have axial 
symmetry are considered. Nomographs 
are given that can ease force calcula- 
tions for the ideal case, and experimental 
test results are discussed. 


Analysis of Electrodynamic Contact 
Forces 


Two methods are used to derive the 
equation for the repulsion force: the first 
is similar to that of Holm, based on a 
direct calculation of the force by consider- 
ing lines of current flow in the constric- 
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tion; the second method is based on 
energy concepts and uses the principle of 
virtual displacement to calculate the 
force. 

In the first method, the contact is 
considered to be a long cylindrical rod with 
a tapered end, making contact on a flat, 
infinite plate. The current flow lines 
are considered as flowing down the axis 
of the rod, through the tapered section and 
contact constriction. Fig. 1 is a cross 
section of this assembly. In the region of 
the tapered section, the flow lines flow 
through a large spherical sector of radius 
A/sinay and terminate on a spherical 
surface of a small spherical sector of 
radius a/sin ap. 


CONTACT ROD 


RADIUS A > 


SPHERICAL SECTOR 
ANGLEQ 5 RADIUS A/SIN a, 


SPHERICAL 
SECTOR 
ANGLE ay 
RADIUS a/SIN ao 


CONSTRICTION 
AREA, RADIUS a 


SPHERICAL SECTOR 
ANGLEa@ 5, RADIUS r 


INFINITE PLATE 


Lc) 
Contract Nonr-1100(8), Office of Naval Research 
Washington, D. C., Setp. 1960. 


+. 
8. DiciTaL AND SAMPLED-DaTa CONTROL Sys. 
tems (book), J. Tou. McGraw-Hill Book 
Company, Inc., 1959. ry 
9. ADAPTIVE SAMPLED Data SySTEMS, B. Wildrow 
Butterworth’s Scientific Publications, London 
England. .. 


In Fig. 1 and the relevant equations th 
following definitions apply: q 


ay =vertex angle of spherical conical se 
measured at vertex 0 from the cente 
axis to the tapered contact surfac 

A =radius of cylindrical rod | 

a=radius of contact construction at poi 
of contact with infinite plate 

6=angle measured about center axis of 
(not shown in Fig. 1) 

y=radial distance of an element of 
ume 7?Xsina dr Xda X dé from 
tex 0 

a=vertex angle of spherical conical sec 
of radius r that defines position 
element of volume 7?XsinaXdrXdd 


x dé 7 
I=total current through spherical conte 
a and contact constriction 
Iz=portion of total current through sphet 
cal conical a 
J, =current density at any point in spherice 
conical a at radius 7 - | 
7=flux density existing at an element c 
volume on edge of spherical @ a! 
radius 7 
u=permeability of contact medium 
ds=surface area of element of volume 2 
radius r perpendicular to radius ling 
r; ds=r?Xsin aXdaXdé 
df =increment of force acting on element c 
volume perpendicular to radius r 


ELEMENT OF VOLUME 
ON EDGE OF SPHERICAL 
SECTOR ANGLE w 
RADIUS r 


CONTACT SURFACE TAPER 


Fig. 1. Cross section of tapered rod making contact on infinite plate ' 
: 
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REPULSION 
FORCE F 


COAXIAL RADIUS M 


CONSTRICTION 
RADIUS a 


CONTACT ROD 
RADIUS A 


CONTACT 
FACE 


CURVATURE C 


1=component of df acting parallel to axis 
of contact rod 

=integral of df; over contact volume equal 
to total electrodynamic repulsion 
force acting parallel to axis of con- 
tact rod 


Uniform J, is assumed at any r from 
rtex 0 of the spherical sector ay. A 
ortion of this current (J;XdS) will pass 
ough the element volume r?Xsin aX 
rXdaXdé and will react electrody- 
amically with the internal B, at this 
oint. The elemental force df produced is 
erpendicular to 7. A component force 
i=sin aXdf exists that reacts parallel 
o the rod axis. The total contact repul- 
ion force F= f° df; Xdv (integral df, over 
olume of spherical ao) reacts on the con- 
act in a direction parallel to the rod 
xis. 
| B,;can be calculated by means of vector 
potential and total circuit current dis- 
However, if current flow 
ough the contact rod is assumed to be 
mimetrical with the rod axis, and the 
eturn path displaced a great distance, 
Ampere’s law can be used by considering 
jhe amount of current flowing through the 
ipherical a. 
Now, B, can be calculated as follows: 
lle pI (1—cos a) 


el (i) 


2ar(sine) 2ar(1—cos ag)sin a 
here 


a=total current X(surface area of conical 
)/(surface area of conical ap) 


i (1— . 
2. (1—cos a) (2) 
1—cos ap 


Now, J; is total current/(surface area 
bf spherical ay at 7): 


Se 
~ 2xr(1—cos ao) 


(3) 


Then, to calculate the repulsion force, 


Fig. 2. Coaxial contact circuit element 


df,=df sin a=B,;XJ,dsXdrXsin a (4) 
where 


df=B,XJdsXdr (5) 


Substituting equations 1 and 3 into equa- 
tion 4 yields 


ul? (1—cos a) 


df, =——_- ———— _ si dr d 
Lf AEN FR cr sin a dr dad@ (6) 
Boe 
in & 
le: i 
= Gee ot 
ih fi An? a ees 
sin % 
ed cos a) sin ad oe 
oo a 
iE do (7) 
f [1—cos ap]? 0 
A 
a8 at 24 (2 a —2 cos ap +1 
Qr = r 2\cos?a9—2 cos ag +1 
pl? A sin ay pl? A 
4 u(S ao a ) Speake (8) 


The repulsion force is found to depend 
on J*, w of the contact material, and on 
the natural logarithm of the ratio of the 
contact radius A and constriction radius 
a. Since the portion of the integral 
that concerns the taper angle ap is unity, 
the force is independent of this taper 
angle. 

A derivation of eaten 8 based on 
energy concepts can demonstrate that the 
repulsion force is independent of an 
arbitrary contact surface curvature. If 
the contact and its return circuit are con- 
sidered as a circuit element, the repulsion 
force can be calculated by considering 
energy changes in this circuit element 
when a virtual displacement of the contact 
is made. 

Fig. 2 is a cross section of the contact 
and current return path formed into a 
coaxial circuit element. The return path 
is shaped as a hollow tube closed at one 
end. The contact is a movable rod placed 


in the center of the tube making contact 
at the closed end. The contact surface 
curvature of the rod is quite arbitrary, 
restricted only to having symmetry with 
the rod axis. The arrows show the direc- 
tion of current flow in Fig. 2, where 


T=total current through contact 

M=radius of current coaxial return path 

A =radius of contact rod 

a=radius of contact constriction, considered 
as a small cylinder 

c=region of contact curvature of an arbi- 
trary shape, but symmetrical with 
the center axis of the rod 


When the contact moves a virtual dis- 
tance Ax, changes of energy in three 
areas occur as follows: 


1. Change of energy stored in the circuit 
element magnetic field identified as W;, 


2. Work done by contact in moving virtual 
distance Ax identified as Wm. 


8. Change of energy supplied by the ex- 
ternal circuit identified as We. 


Then, by the eS of energy con- 
servation, 


We- Wm=W: (9) 


The contact circuit element is driven 
by a current source, since the impedance 
of a contact is very low compared with that 
of the rest of the circuit. This means that 
changes in contact impedance cause little 
change in current; hence the current 
through the contact can be assumed con- 
stant during the virtual movement Ax. 

If the virtual movement Ax occurs in 
time AT, the electric energy into the con- 
tact element during this time is 


We=S LE d (10) 


where 


E=induced voltage across contact element 
caused by changing flux linkages 
when virtual movement is made 

J =instantaneous line current through con- 
tact assumed to be constant over 
time Af 


Then, 


: 
(LI) fs it 

Wr= I d—dt= I*dL LIdI 
: IPs (di) ie 


Since J; =J2.=a constant, 


= J?(L2—L1) =I(A2—1) (11) 


where 


Ii, \:=inductance or flux linkages of the 
contact circuit element at zero refer- 
ence time 

=inductance or flux linkages of the 


contact circuit element at time AT 


Lx, d2 


The change of internal energy W; of the 
circuit element is proportional to the 
difference in its inductance before and 
after the virtual contact movement. 
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NOMOGRAPHS FOR THE CALCULATION OF 4 
CONTACT ELECTRODYNAMIC REPULSION : 


FORCES 


SS 
ou 
Ore) 
x< 
fo) 
(ey) 


Noo 
folerey 


GRAD 
O000 


12.0 


Fig. 3. Repulsion 
force nomographs for 
currents, in amperes, 
up to: A—3,000; 
B—10,000; cC— 
25,000 
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I = {INSTANTANEOUS CURRENT THROUGH CONTACT "a SPOT IN AMPERES 


5.0 


cuit element caused by virtual contact 
displacement during time AT is 


20 


(B) 


fs B12 D 
F=36x 1081? Ln( avian x 7) 


H 
SCALE FACTOR = peers) 


100,000 


ft 
| 
: 
, 
Ax =virtual movement of contact duri ig | 
time T 


Substituting equations 11, 12, and 138} 
into equation 9 yields 


1 1 


(14)) 


lating the flux linkages of the contact co- 
axial circuit element before and after a 
virtual movement Ax. Since the system is 
coaxial, the simple equation commonly 
applied to calculate total flux linkages per’ 
‘unit length of a coaxial wire can be used. 


Ne af oe ie 
ae 4] (1 


where M is the inside radius of the coaxial | 
return path and + the radius of the: 
inner conductor. The term pl/2rX1/4 
defines the internal flux linkage of the 
inner conductor. This internal flux link- 
age is independent of conductor radius if 
a uniform current density over a cros 


5 
The problem reduces to one of calcu-; 
| 
) 
I 


we 


1 f ih owen section is assumed.? 
ae T(x) (12) Wm=Jo F dx=FAx (13) When the contact of Fig. 2 is moved 
where distance Ax, the flux linkages are changed 
The mechanical energy out of the cir- F=repulsion force in two areas: the portion of flux linkage 
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rounding contact rod radius are de- 
mased by an amount plAx/2r[(In M/- 
+1/4]. The flux linkages of the 
lindrical conductor comprising the 
tact constriction are increased by an 


He flux linkages of region C, defining the 
Wbitrary face curvature of a symmetrical 
fntact, remain unchanged during a 
fovement of Aw. 

The net change in flux linkages with 
is movement are 


M 1 M 
| (in 2) 44—(in t)=i1 (16) 


um ( A x) plAx. A 
in = = In 


pitas. A 
(A2—A1) = In — = FAx (17) 
alt 4 
4a a 


quation 17 is the same as equation 8 de- 
ed by the first method. 


ractical Calculations 


} Equations 8 and 17 demonstrate that 
ine repulsion fjorce is dependent on the 
pntact constr ction radius a. Holm! 
ows that the magnitude of a can be 
pproximated by considering that the 
Iressure between contacts causes the 
idging constriction to be in a state of 
lastic deformation that supports the 
bad. It is assumed that the contacts 
puch in only one point and that the area 
If contact is circular. The equation re- 
Liine this contact area to a material con- 
ant is given by Holm‘ and Shobert® to 


(18) 


a=constriction area, inches? 

7,=contact material hardness, pounds per 
inch? (see references 4 and 5) 

1=contact pressure in pounds 


f the area A, is assumed to be circular, 
hen the radius a of the constriction is 


(19) 


The nomograph calculation methods 
hat follow use a constriction radius cal- 
ated from equation 19 and are based on 
equations 8 or 17, Use of these equations 
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Fig. 4. Calculated comparison of influences on repulsion force as functions of current. 


A—Influence of H, for contacts of 1-inch diameter. 


B—Influence of contact diameter; H = 


100,000 psi 


requires that the magnitude of constric- 
tion radius be determined by the initial 
contact pressure, and that the radius re- 
mains at this magnitude until the contacts 
separate or “‘blow apart.” Actually, this 
idealized condition is not realized in prac- 
tice. However, experiment shows that 
the idealized momograph calculations can 
be modified to give meaningful results. 

If practical measurement units and the 
magnitude of a from equation 19 are sub- 
stituted into equation 8 or equation 17, 
the equation from which the nomographs 
in Figs. 3(A), (B), and (C) were designed 
is obtained. This equation is 
2\/ nH. D 

vos 5) (20) 


F=5=30%10-474In( 


where 


F is expressed in ounces 

S=applied contact pressure ounces 

P=net contact pressure equal to difference 
between S and F, ounces; P=S—F 

H is expressed in psi (pounds per square 
inch)!» 

D=contact diameter, inches : 
T=instantaneous current through contact 
constriction or a spot, amperes 

The nomographs were designed on the 
basis that the contacts separate when 


F=S. 


Use of Nomographs 


The nomographs given in Figs. 3 (A), 
(B), and (C) are designed for various con- 
tact diameters of materials with an H of 
100,000 psi. The force is determined by 
laying a straight edge from the value of 
contact diameter to the value of instan- 
taneous current through the contact con- 
striction or aspot. The force required to 
keep the contacts from separating is read 
on the centers cale. 


It should be noted there is a separate 
calibration called a scale factor, which 
is used for making calculations involving 
materials with an H different from 100,000 
psi. The determination of H for mate- 
rials can be made from superficial Rock- 
well hardness readings as shown by Sho- 


bert.5 The scale factor (SF) is calculated 
by 

| H 
SP Sa (7) oY 


This magnitude is located on the SF 
scale on the appropriate nomograph. 
The dimension is transferred with dividers 
to the D scale. The SF dimension is 
added or subtracted to or from the loca- 
tion of the required D, This new loca- 
tion is used to make calculations for the 
required D. 

As an example, suppose it is desired 
to determine the force on a copper con- 
tact with an H of 175,000 psi, 1-inch 
diameter, carrying an instantaneous cur- 


rent of 5,000 amperes. The SF calculated 
by equation 21 is 
1. (175,000) 
B= =) .278 22 
S= 3)" Gonacoy @ 


With the use of a pair of dividers the SF 
dimension +0.278 is added to the 1-inch 
location on the D scale, which is the point 
to use to calculate the contact repulsion 
force for a material of an H=175,000 psi, 
one,inch in diameter. A line is drawn 
from this point to the instantaneous cur- 
rent of 5,000 amperes, after which 
50 ounces is read on the force scale. 

Fig. 4(A) shows a family of curves giv- 
ing calculated contact repulsion force as a 
function of current for various values of 
H. The curve shows that a variation of 
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Table |. Experimental Results 
Force,* Ounces 
Current, Per Cent 
Amperes Calculated Measured Difference 


5,500 65 SO erry o. 


* Contact pressure required to keep contacts from 
separating when carrying given current. 


H of 4 to 1 causes a change in repulsion 
force on the order of 20%. 

Fig. 4(B) shows calculated repulsion 
force as a function of current for contact 
diameters of 1/8 inch and 1 inch of a 
material with an H of 100,000 psi. 
Examination of the curves shows that this 
variation in diameter causes a repulsion 
force change on the order of 60%. 


Experimental Results 


Table I lists the results of tests per- 
formed on 1-inch-diameter copper con- 
tacts. It is evident that the contact 
pressure needed to hold the contacts closed 
varies by about 20%. Part of this varia- 
tion can beattributed toerrors in measure- 
ment of constriction radius a. If contact 
conditions are ideal, so that the constric- 
tion area is circular and in a state of 
plastic deformation, measurement of con- 
tact pressure should give a reasonably 
accurate value of constriction radius. 


However, as Holm‘ shows, the constric- 
tion region may be partially in a state of 
elastic deformation as well as being plas- 
tically deformed. This deviation from 
the ideal case in a practical contact causes 
errors in calculating the constriction 
radius a by equation 19. This error would 
be such that a smaller value of repulsion 
force is calculated by the nomographs 
than actually exist. 

The magnitude of contact mass and 
pressure-spring rate can also cause a 
variation of measure repulsion force from 
thatcalculated bythe nomographs. When 
the repulsion force is exerted on the con- 
tact, the contact will start to move in ac- 
cordance with the dynamic properties of 
contact massand pressure-spring constant. 
The movement will reduce the size of the 
constriction radius a, which will cause a 
greater repulsion force at a given current 
than wouldnormally exist. This increased 
force will cause an increase in contact 
acceleration and movement, leading again 
to an increase in repulsion force. The net 
result wili be to part the contacts at a 
lower value of current than is calculated 
by the nomographs. 

Even though these factors cause force 
magnitudes calculated by the nomo- 
graphs to disagree with experimentally 
measured values, the nomographs provide 
a good tool to estimate contact pressure 
force requirements. Experience with a 
variety of contact structures demonstrates 
that for any given class of contact struc- 
ture a multiplying correction factor can 
be applied to the results. For example, a 
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B. ROLSMA 


NONMEMBER AIEE 


HIS PAPER is a discussion of the 

principles of steady-state heat transfer 
as they apply to the field of electric heat- 
ing and, as such, no solutions are derived, 
but rather the governing differential 
equations are presented and references 
containing the solutions for various bound- 
ary conditions are cited. In most of the 
references given, the work was concerned 
with the temperature distribution in elec- 
tric coils or other electric apparatus where- 
in the generation of heat was an un- 
wanted phenomenon and not, as is consid- 
ered here, the prime function of the de- 


tions to basic differential equations, they 
are directly applicable to any other body 
for which the basic equation and the 
boundary conditions hold. 

An attempt has been made to divide 
the discussion into essentially two parts; 
the first part deals with the heat transfer 
within electric heating elements them- 
selves, and the second part is devoted to 
heat transfer from the surfaces of the 
elements. Since the temperature dis- 
tribution within an electric heating ele- 
ment depends upon the heat transfer 
rate from its surfaces, the two are not 


t <a 


multiplying factor of 1.2 applied to the 
calculated data shown in Table I woule 
insure that the calculated contact p 

sure determined would keep the “- 
closed. ‘ 


Conclusions 


4 
The analysis has provided theoretiea 
insight into the nature of the contact 
repulsion force to the extent that this foree 
is independent of contact face curvatt 
with contacts shaped symmetrically about 
the axis. In addition, nomograph cal- 
culation aids have been provided an 
estimate repulsion force magnitude. H 
ever, more study is needed to account for 
the variance between measured and calecu- 
lated force magnitudes. The foremos 
problem is that of analyzing in detail the 
effect of contact mass and pressure-spring 
constant on the repulsion force. 
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engineer applying purchased electri 
heating elements is interested in dissipat- 
ing the developed energy from the s 
faces of the elements and not in the inter- 
nal heat transfer problems. 


Exact Solutions for Heat Transfer 
in Electric Heating 


$ 
q 
? 
‘ 
The differential equations for the 
steady-state temperature distribution if 
an electric heating element are deriv 
by performing a heat balance upon — 
{ 
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fferential volume in the solid material. 
WSsuming that the medium is homo- 
#neous and neglecting any changes in 
iS thermal conductivity due to tempera- 
€ variations within the heating element, 


he heat balance on the differential volume 
telds 


O7t O% d3% 
sit pty) t2""=0 (1) 


If the internal heat generation per 
nit volume, the q’”’ term, is considered 
> be constant, this equation is exactly 
oisson’s equation. 
Solutions to equation 1 for simple geo- 
netric configurations with fixed surface 
smperatures and constant g’’’ are well 
mown. Rogowski! has shown that the 
perature distribution within infinitely 
ride plates, infinitely long cylinders, and 
pheres has a parabolic form. Jakob? 
eated the infinitely long rectangular bar 
ith these same boundary conditions. 
he results of his work expressed as a 
atio of the temperature distribution in 
n infinitely wide plate are summarized 
n reference 3. 
The assumption of a negligible change 
n the thermal conductivity with tempera- 
re as used in the derivation of equation 1 
oupled with the use of a constant g’’’ 
etm in producing the solutions cited is 
uite realistic for electric resistance 
eating element materials. The Wiede- 
mann-Franz-Lorenz equation indicates 
hat thermal conductivity and electric 
onductivity are related by a simple 
xpression of temperature multiplied by 
Lorenz number. Since the Lorenz num- 
per does not exhibit large changes with 
jemperature variations, thermal conduc- 
ivity and electric conductivity are essen- 
lially related by an expression of a con- 
itant multiplied by temperature and, 
hence, the ratio, g’’’/k, in equation 1 may 
be expected to follow this same relation- 
thip. The problem then is to select a 
value of g'’’/k in the range of tempera- 
ure encountered in the problem. 
Although it is recognized that both 
hermal conductivity and heat generation 
n electric resistance heating elements vary 
with temperature, solutions to equation 1 
nave been developed considering that the 
"’ term alone varies with temperature. 
akob‘ generated solutions to equation 1 
or the infinitely wide plate, the infinitely 
ong cylinder, and the sphere with fixed 
surface temperatures and a linear varia- 
ion of q’’’ with temperature, while 
iggins® investigated the temperature 
listribution in an infinitely long rec- 
angular bar under these same condi- 
ions. 
When the changes in thermal conduc- 
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tivity with temperature variation cannot 

be negelected in a homogeneous material, 

the steady-state heat balance on the 

differential volume yields the relationship 
E o% 4 


Ok] ( ot \? (de \2 ot\? 
a tte) +e) row 

This equation is clearly nonlinear. 
However, the use of the transform pro- 
cedure introduced by Kirchhoff® in 1894 
will reduce equation 2 to a linear equa- 
tion similar in form to equation 1. Ap- 
plications of this procedure are to be found 
in the discussion by Jakob of Higgins’ 
paper’ and in reference 8. 

The transform method is applicable 
to the cases of varying thermal conduc- 
tivity with a constant g’’’ term or with a 
temperature-dependent g’’’ term. When 
it is considered that both thermal con- 
ductivity and gq’ vary with tempera- 
ture, both & and q’’’ are expressed in 
terms of the transform variable and sub- 
stituted in equation 2. This substitution 
will reduce equation 2 to Poisson’s equa- 
tion. On the other hand, when thermal 
conductivity alone varies, the solutions 
obtained for equation 1 with a constant 
q’’’ term may be directly applied to 
equation 2 by introducing the transform 
variable into these solutions provided 
the boundary conditions are those of 
fixed surface temperature or thermal 
conductivity multiplied by the tempera- 
ture gradient at the surface. 

The heat balance for the more general 
case of heat transfer with thermal 
conductivity varying with both tempera- 
ture and direction in the differential 
volume yields 


Q MMO MOEN \ dif ads 
—( ke — )4+—( ky — )4+=( 2 = 
oie =) +2( ‘i = +2( =) if 


. a"=0 (3) 


Reference 8 contains a solution for this 
case considering that the g’’’ term varies 
linearly with temperature and, further, 
that the heat transfer rates at the four 
surfaces of the rectangular bar examined 
vary, i.e., both the surface heat transfer 
coefficients and the sink temperatures 
are different for the four surfaces. 

The exact analytical solution of the 
general case of heat transfer described 
by equation 3 with variable heat genera- 
tion and the boundary conditions of heat 
transfer with different surface coefficients 
and sink temperatures on all faces of the 
body for even the more simple geometric 
shapes is a very difficult and time-con- 
suming task. Quite often it is not even 
possible to obtain the solution by exact 


analytical means and finite difference 
techniques must be used. However, be- 
fore leaving the area of exact analytical 
solutions, the work of Carslaw and Jaeger® 
must be mentioned for it is the most 
complete compilation of mathematically 
rigorous solutions to the problems being 
considered here. 


Approximate Solutions for Heat 
Transfer in Electric Heating 


In the numerical analysis of heat 
transfer problems, the differential equa- 
tions 1, 2, and 3 are replaced by differ- 
ence equations of the general form 


E tn4 Az,y,2) — 2l(2,y,2) HE (oe Ary 2) 4 
ri (Ax)? 


pean see) Hewson | 
(Ay) 

ke L(x yet Az) — 2k (we y,2) th (2,y,2- Az) aa 
(Az)? 


q’’'(@,v,z))=0 (4) 


Schneider” has applied numerical tech- 
niques to yield solutions for the tempera- 
ture fields in electrical coils, while Dusin- 
berre! presents a very complete case for 
the solution of heat transfer problems 
in general by numerical analysis. 

While the relaxation method originated 
by Southwell? and initially applied to 
problems in the field of heat transfer by 
Emmons!’ is somewhat tedious, it will pro- 
duce acceptable results provided the 
body subdivisions, Ax, Ay, and Az, are 
small. As Ax, Ay, and Az become ex- 
tremely small, the results of the relaxation 
process approaches the exact solution for 
a given problem, and, in the limit, as 
Ax, Ay, and Az approach zero, equation 4 
becomes equation 3 exactly. However, 
as Ax, Ay, and Az become smaller in. a 
given region, more and more finite differ- 
ence equations are introduced. Some 
compromise between the desired accuracy 
of results and the time or labor involved 
in obtaining a solution is then in order. 

When a given problem is first encoun- 
tered, there is no sure way of judging the 
errors introduced by selecting a certain 
body subdivision. The only course open 
is to select one subdivision size and to 
obtain an initial solution. The analysis 
is then repeated with a subdivision size, 
say one-half the previous one. A wide 
difference in the results of these two solu- 
tions would indicate that a further re- 
duction in subdivision size is required; 
while a very small difference in results 
indicates that a larger subdivision could 
possibly be used. It is only through ex- 
perience of this nature that a basis for 
selecting subdivision size can be derived. 
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The relaxation method does, however, 
have one significant drawback in that it 
is not readily adaptable to automatic 
sequence-controlled digital computers. 
This disadvantage is not a serious con- 
sideration in organizations without digital 


computing facilities, but where such 
equipment is available, other methods are 
being employed to solve the finite differ- 
ence equations. 

One of the methods frequently used in 
conjunction with digital computers to 
solve finite difference heat flow equations 
is the iteration method.'4 Although the 
iteration process, when performed man- 
ually, is generally much slower than the 
relaxation process, general programs 
may be written to direct a digital com- 
puter to automatically and rapidly per- 
form the calculations requ'red by the 
iteration process and free the analyst 
of this task. 

The speed of the digital computer and 
its ability to store and process much data 
somewhat eliminates the need for the 
compromise between desired accuracy of 
results and subdivision size mentioned 
before. Using the digital computer, it is 
possible to evaluate quickly the effect of 
subdivision size on the results and to 
actually “home in” on a solution by re- 
peatedly reducing the size of the sub- 
divisions and rerunning the problem 
through the computer. The only appar- 
ent restriction would be computer memory 
size and the time required to complete 
the data input forms necessary for the 
program. 

The digital computer is then a power- 
ful tool, but a word of caution: The re- 
sults of the calculations can be no better 
than the data used. Again, in the case of 
uncertain data, the speed of the digital 
computer does offer an advantage in that 
a problem may be rapidly solved with first 
the high value of the uncertain data and 
then with the low value of the data. A 


30 


Fig. 1 (left). 
Various sizes of 
General Electric 
Company  cart- 
ridge - type 
heaters 


Fig. 2 (right). 
1,000-watt water 
immersion heater 


comparison of the two solutions would 
indicate whether a better value for the 
data is required. 

In the derivation of a number of the 
previously cited solutions, it was assumed 
that the surface temperatures of the 
heating elements were fixed and invariant 
over the individual surfaces. This may 
not always be the case in practice, for the 
surfaces of a heating element and even dis- 
crete portions of these surfaces will 
assume whatever temperature may be 
required to dissipate the heat generated. 
If the heat generated within the element 
cannot be dissipated at temperatures be- 
low the melting point of the material then 
the element will fail. Thus, an element 
designed for a certain ambient condition 
and installed under a different condition 
may fail very quickly in service if the in- 
stallation conditions have an adverse 
effect on the dissipation of heat from the 
surfaces of the element. It is imperative, 
then, that the design engineer be informed 
of all the conditions of installation when 
custom-designing heating elements, and 
that the equipment design engineer, 
when applying ‘‘off-the-shelf’’ heating 
elements, be fully aware of the design 
conditions of these heating elements. 


Heat Transfer from Electric Heating 
Elements 


All of the various modes of heat trans- 
fer, conduction, free and forced convection, 
and radiation, may be active to one extent 
or another in removing heat from the sur- 
faces of electric heating elements and in 
some cases all may be important. As an 
example of this, consider the tubular-type 
heating element found in the modern 
electric range. This unit is fabricated 
by centering a helically wound section of 
resistance wire within the element sheath 
and then filling the space remaining in- 
side the tube with an electrical insulating 
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powder such as magnesium oxide. 
nally, the tubes are swaged or rolled tor 
duce their diameter and to compress 


powder. The heat generated within a 
resistance wire is transferred to tt 
outer tubing by conduction through fi 
insulating powder and, in addition, f 
radiation between particles. Heat 
dissipated, in turn, from the outer tubint 
or sheath, by radiation and free conte 
tion. 

In some cases, conduction is the pri 
cipal heat transfer mechanism removiti 
heat from the outer sheath of heati 
elements. Low-temperature, 750 F é 
grees Fahrenheit) maximum, cartrid 
heaters of the type shown in Fig. 1 aa 
inserted in holes in die blocks, pre 
platens, molds, and other bodies whid 
must be heated. These heaters are a: 
signed to fit snugly in a hole of a specifi 
size and to dissipate their energ 
through the gas film in the clearance sp 
between the heater and the sides of : 
hole to the body being heated. Sina 
radiation is not very effective at these lox 
temperatures, almost all of the heat wi 
be transferred by conduction through th! 
gas film. If by some chance the size « 
the hole is increased slightly so that th 
clearance is twice the design value, th 
temperature differential between tt 
heater sheath and the surfaces of the ho! 
will be approximately twice the desig 
This increase in temperatum 
differential may cause the temperatum 
within the heater to rise beyond th 
maximum permissible. 

Immersion-type heaters, Fig. 2, i 
designed to dissipate their heat Bb 
natural convection to liquids and, heneé 
any installation condition which inhibit 
the free circulation of the liquids over tH 
heater element surfaces is to be avoide¢ 
In regard to this, locations which plac 
immersion-type heaters in close proximit 
to the sides of the containing vessel ¢ 
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par the free surface of the liquid being 
pated are poor. 

It should also be recalled that the liquid 
hysical properties—viscosity, thermal 
bnductivity, specific heat, and density— 
e important factors in determining the 
irface coefficient of heat transfer for both 
reed and free convection. Thus, an 
nmersion heater designed for heating 
ater to a certain temperature may not 
e suitable for heating other liquids to this 
ume temperature, since 4-to-1 variations 
the natural convection coefficient be- 
ween different liquids is not uncom- 
on. 

Attention must be directed to the 
inciples of forced-convection heat trans- 
+ and fluid flow when heaters of the types 
own in Figs. 3 and 4 are used to heat 
ases flowing in ducts. When heaters 
e to be installed in ducts or boxes 
nrough which gases flow, configurations 
ilar to those encountered in the de- 
gn of heat exchangers are involved. 
uch information on the flow character- 
tics and heat transfer rates from bare 
md finned tubes is to be found in the 
echnical literature. 

The flow information obtained from 
ests of heat exchanger configurations 
an be directly applied whenever these 
me configurations occur in the applica- 
on of electric heating elements. Thus, 
he bypass effects found in heat exchang- 
ts will also be present in electric heating 
pplications. The spacings required to 
btain good mixing between rows of 
heat exchanger tubes will also be required 
between rows of electric heating elements 
nd the pressure drop through rows of 
lectric heating elements may be cal- 
ated from correlations for the pressure 
lrop in heat exchangers. 


Marcu 1961 


On the other hand, the heat transfer 
data derived from heat exchanger studies 
cannot be directly used for determin- 
ing electric heating element designs. 
Schmidt and Wenner!!8 have shown that 
the local surface heat transfer coefficients 
vary along the circumference of a 
cylinder when flow is normal to its axis. 
This inherent variation in the local co- 
efficient around the surface of a cylinder 
can readily, and indeed does, create 
localized hot spots on the surface of 
cylindrical heating elements in cross- 
flow. This variation in coefficient is 
neglected in correlating data from heat 
exchanger investigations and the values 
of heat transfer coefficient reported are 
average values taken over the entire sur- 
face of the tubes. These coefficients are 
then not suitable for electric heating 
element design, for it is the coefficient 
based on hot-spot temperature which 
must be used to determine the electric 
heater rating. 

One additional consideration, not pres- 
ent in the design of heat exchangers, does 
exist in selecting and arranging electric 
heating elements for heating gases, or 
for that matter, any fluid flowing in a 
duct or pipe. When large temperature 
rises are required in the fluid flowing 
through the heater bank, it may be neces- 
sary to derate or reduce the total power 
input to the heaters in the last few rows 
of the heater bank so that their tempera- 
tures do not exceed the maximum allow- 
able. 

When temperatures in excess of 700 F 
are encountered in electric heating, heat 
transfer by radiation begins to play an 
important role. While it is well recog- 
nized that the temperature level of opera- 
tion is a prime factor in establishing the 


rate of heat transfer by radiation, there 
are other factors which must also be con- 
sidered. The geometry of the heating 
element, the geometry of the heat re- 
ceiver, and their positions relative to each 
other and to any other surfaces viewed by 
either the source or the receiver are fac- 
tors which also influence the rate of 
radiation heat transfer. In addition to 
these geometric considerations, the surface 
characteristic, emissivity, must be ex- 
amined for not only the energy source 
and the ultimate receiver of this energy, 
but for all surfaces in the radiation sys- 
tem. 

Much information is available con- 
cerning the geometry factors in radiation 
heat transfer.!7—!9 Emissivity values for 
both unoxidized and oxidized metallic 
surfaces and for nonmetallic surfaces are 
to be found in references 20 and 21. 

The solution to problems involving 
radiation heat transfer can become quite 
complex when more than three surfaces 
are involved in the interchange. Not 
only will the consideration of additional 
surfaces make the problem more complex, 
but the simplest 2-surface problem be- 
comes complex when variation in the 
temperatures of these surfaces makes it 
necessary to subdivide them into zones of 
uniform, or substantially uniform, tem- 
perature. Oppenheim’? presents a net- 
work analysis method for solving the 
diffuse radiation problem within an en- 
closure formed by grey surfaces. Hottel”? 
presents a slightly different approach to 
this same problem. These methods, for. 
the steady-state case, result in a series 
of linear simultaneous equations which 
can be solved either with a digital com- 
puter, if one is available, or with a desk 
calculator, for up to eight equations, 


Fig. 3 (left). Edge-wound resistor heating unit for use in forced- 
convection heating 


Fig. 4 (above). 


Finned tubular heaters 
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Fig. 5. Electric resistor heating element for furnace service 


using Crouts?® method as suggested by 
Hottel. 

Heating elements for high-temperature 
service in ovens and furnaces are made 
from both metallic and nonmetallic elec- 
tric resistance materials. Nonmetallic 
elements are primarily straight cylindrical 
shapes with low-resistance sections at the 
ends where the elements pass through 
the furnace or oven walls, while metallic 
heating elements are frequently cast or 
formed in sinuous shapes. Figs, 5 and 6 
show metallic elements formed using 
resistance alloy ribbon and molybdenum 
rod respectively. Also shown in Fig. 5 
are the low-resistance terminal ends to be 
passed through the furnace wall. 

When designing the metallic resistors, 
the spacing between adjacent loops must 
be made large enough to permit the heat 
generated within the material to radi- 
ate freely out into the furnace. If the 
spacing is too small, the ribbon tempera- 
ture required to dissipate the energy 
developed will become unreasonable. 
The length of each loop is controlled by 
the ability of the material to support 
itself at high temperature. If the loops 
are too long, sagging or excessive elonga- 
tion of the elements will result. 

Attention to the details of supporting 
the resistance heating elements in the 
furnace may reduce the operating tem- 
perature of the element considerably. 
When ceramic supports are employed, 
localized hot spots may develop because of 
the radiation screening or shading effect of 
the support and the inability of the 
ceramic material to conduct heat away 
from these local hot spots. Metallic 
hooks are preferred for supporting the 
elements, since they are usually much 
smaller than ceramic supports and, hence, 
do not screen the view from the element to 
the furnace as much, and have, in addition, 
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the ability to conduct the heat away from 
a hot spot more readily then ceramic 
materials. 

With adequate support, the only limita- 
tion on the temperatures attained by 
electric heating elements is the melting 
point of the element material. However, 
since the oxidation rates at these elevated 
temperatures would normally result in a 
very short element life, protective or inert 
atmospheres are introduced into the 
furnace, or the furnaces are designed for 
operation with a vacuum atmosphere. At 
first glance, the requirement of a protec- 
tive or inert atmosphere for reasonable 
element life would appear to be a deter- 
rent to the use of electric heat, but it 
must be remembered that many of the 
processes which require the higher tem- 
perature. levels also require that the 
processing be done in a_ protective 
atmosphere. Actually then, electric heat 
has an advantage over other methods 
of heating in this case, since electric 
elements can operate in protective atmos- 
pheres without contaminating the atmos- 
phere, while other methods of heating 
require the use of an envelope or muffle to 
separate the heat source from the protec- 
tive atmosphere and the material being 
processed. Further, whenever a muffle 
or envelope is placed between the source 
of heat and the eventual receiver of heat, 
the net rate of radiant interchange at 
a specified temperature level is reduced. 


Conclusions 


In summary, the heat transfer problems 
of the electric apparatus designer and 
the electric heating element designer are 
closely related, for both are concerned 
with the same basic differential equations 
within the heat generating regions of their 
respective equipments and with removing 


Fig. 6. Molybdenum resistor heating element for high-temperature 


Sea 


furnaces 


this heat to the equipment surroundings 
with as small a temperature different 
as possible. The only significant diffe 
ence between the two is that the elec 
apparatus designer considers the Joule 
heat as an unwanted quantity and se 
to reduce it to an absolute minim 


signer is constantly striving to incre 
the heat generated to the maximu 
possible without sacrificing element life. 


Nomenclature 


k=thermal conductivity; kz, Ry, 
thermal conductivity in x, 
direction 

q''’=internal heat generation per unit vol- 
ume 

t=temperature 

x, Y, 2=co-ordinate axis directions 7 

Ax, Ay, Az=increments in the x, y, z 
directions ; 
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solution; then, to show how a computer 
can aid standard servomechanism pro- 
cedures to effect an optimum design, with 
the tedious portion of the work left to the 
computer. 


System Type 


The system described, which is char- 
acteristic of this type, saturates at a 
relatively small error signal. Thus, the 
discussion is confined to system perform- 
ance with small signals. Saturation is 
effected at small signal levels since, to 
minimize positioning time, proportional 
action is desired only for the minimum 
distance from the set point necessary to 
achieve the required positioning accuracy 
and performance. Thus, the control sys- 
tem designer is concerned with optimizing 
the performance during that portion of the 
positioning cycle when the motion ap- 
proaches its final position. The time that 
is required to position after the controller 
leaves its saturated state is called its 
settling time. 


How to Determine Equations 


Essentially, the method for determining 
the equations of motion makes use of the 
system’s frequency response to effect an 
analytic expression for the transfer func- 
tion. The plotted data from a frequency 
response test is the well-known Bode plot 
where both response of the system in 
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decibels and phase of the system in 
degrees, compared with input, are plotted 
against the logarithm of frequency. The 
plot usually covers three to five decades 
on a setilog graph. 


USE OF FREQUENCY AND TRANSIENT 
RESPONSE 


The frequency response may be used 
to indicate the degree of stability of a 
given system and give an indication of 
compensating networks which will im- 
prove the performance. An analytic ex- 
pression of this plot, while not necessary 
for determining stability, is an aid to 
stability improvement, and it enables the 
designer accurately to pinpoint the type 
of response the system will have to any ~ 
input whatsoever. 

The transfer function might also be 
obtained from an analysis of the transient 
response of the system, but this was found 
by trial to be less accurate than the use of 
the frequency response. This conclusion 
is supported by the work of others.! 


Form of Transfer Function 


The usual form of the transfer function 
is a fraction consisting of a group of first- 
and/or second-order factors in both the 
numerator and the denominator. From 
the slopes of the attenuation plot and the 
corresponding phase angles, a person with 
only a little practice can approximate the 
total order of both the numerator and 
denominator of the sought-after transfer 
function. Even if the exact order of a 
transfer function cannot be estimated 
accurately the first time, by using the 
method described herein, the order of 
the transfer function becomes apparent as 
the steps in this method are progressed. 


Program Steps 


Steps are given herewith for determin- 
ing an analytic expression of a Bode plot 
obtained from test data: 
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i: 0 START 
100, 
| —30 
READ IN: Wi,W2... 
BOF is PARAMETER BEING 
Db: | 
60 —— 's. 
@ LAG | _lg0 ABOVE WS 
AMP(Db)40 DLAG 
| = 120 
20 - ; 
| 150 
| 
‘ | 
= = — ——_——|80 
“kA ae ; 2I0 OPERATE ON THE NEXT 
Ol eke 8 10 100 1000 PARAMETER AND ITS 


W (RAD./SEC) 


1. By observation, assume the Laplace 
transform of the transfer function of the 
Bode plot to be of the form 


KG 
(s/tu+1).. .(s/tint+1)(s?/on? + 
, 2pxS/wut+1).. .(5?/e2n? +2 pons /won+1) 


~ (8/4+1).. (8/t3n-+1)(s2/eou?+ 
2ps1S/e4y+ 1) 6 oie (s?/w4n? +2 psns/wan+ 1) 


2. From the shape of the plot, estimate the 
Taps, Which would be the break frequencies 
of the straight line approximations of the 
Bode plots and estimate the weg’s and the 
pagp’s which are the natural frequencies 
and damping factors, respectively, of each 
second-order factor. To simplify the com- 
puter program, estimate these values for 
the first approximation as slightly lower 
than the assumed value. 


3. Program the computer to calculate 
the values of KG for the assumed values of 
the parameters 728, wag, and peg at various 
discrete values of w, and compare the 
calculated values with the observed data. 
Also determine by observation which 
parameter (i.¢., tag OF wags OF pag) has 
the most weight at each discrete w chosen. 


4. Set a reasonable limit on the maximum 
allowable deviation between the observed 
data points and the calculated values : 
and, if the error is outside of the reasonable 
limit, program the computer to increment 
to a larger value the magnitude of the rag 
Of Wag OF pag Which has the most weight 
at the point where the excessive devia- 
tion occurs. : 


5. By an iterative process, repeat the 
preceding two steps, using the incremented 
parameters until the deviations between 
the calculated data points and the empirical 
data points are within the specified toler- 
ance. If the magnitude of deviation is 
outside of the tolerance and is increasing, 
the computer must be programmed to stop 
the iterative process and increment the 
parameter which caused an increased 
deviation to the next lower value. 


6. Program the computer to store the 
best values of parameters as determined 
by step 5, and to use these parameters to 
recalculate and print out the values of 
KG selected in step 3. 


7. Plot the transfer function of the 
computer-calculated analytic expression of 
the transfer function, and compare it with 
the plot of the empirical data. By ob- 
servation, the engineer will be able to 
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Fig. 1 (above). Bode 
plot of uncompen- 
sated system 


Fig. 2 (right). Com- 
puter flow diagram 


determine what changes, if any, should be 
made in the assumed form of the analytic 
expression, and steps 1 through 7 can be 
repeated until the data from the calculated 
analytic expression and the empirical data 
ate within the tolerance specified. 


These steps proceed in a manner more 
tedious than one would care to try without 
the aid of an electronic digital computer. 
This system, however, is an easily under- 
stood, staightforward, and speedy method 
of obtaining the transfer function equa- 
tion, or an analytic expression from an 
empirical data plot. 

An interesting side comment which can 
be made is that, before this particular 
method was devised, several attempts 
where made to adapt six different stand- 
ard curve-fitting procedures to the set of 
data obtained from the frequency response 
test. These attempts were performed by 
mathematicians, including a representa- 
tive of a well-known computer organiza- 


.W;, ASSOCIATED WITH THE 


EMPIRICAL DATA OF 20 LOG KG FOR THE 


CALCULATE 20L0G KG USING THE ASSUMED 
VALUES OF THE PARAMETERS AND DETERMINE 
THE DIFFERENCES BETWEEN THE 

CALCULATED AND EMPIRICAL DATA. 


CHECKED; AND THE 


ARE THE EMPIRICAL DATA 
MINUS THE CALCULATED 


DATA WITHIN THE MAXIMUM 


ALLOWABLE TOLERANCE 


RESTORE 
PARAMETER 


tion, who were well-versed in compute 
techniques; but who were not familia 
with servomechanism design procedures 
The method herein described, which w 
finally used, was proposed and set up b 
the author, an engineer fairly well 
versed in servomechanism design te 

niques, but only vaguely familiar wit 
computer technology. 


This example serves to point out 
fuller use of a digital computer can 


60 


90 
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Fig. 3. Bode plot =20 
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. 4. Transient response of compensated 
system 


Pen in practice to a great extent and 


r a long time in the analog computer 
ald. 


lotting Attenuation Data 


The seven steps previously outlined 
ere applied to the data taken from a 
Dsitioning table used on a layout drilling 
achine. The attenuation data are 
nown plotted in Fig. 1, in a standard 
anner on 5-cycle semilog graph paper. 
includes the transfer function of a d-c 
otor, driven by a thyratron controller in 

open-loop configuration without any 
bmpensating networks in the loop. The 
sst data and the data calculated from the 
bmputer program are a single curve in 
ig. 1 because the difference between the 
Wo curves was only a fraction of a decibel 
any point. 


esults Obtained 


A flow diagram of the computer pro- 
fram is shown in Fig. 2. The results 
fescribed were obtained after four com- 
futer runs. The time spent by the en- 
ineer in setting up the data input for 
fach run was never more than 7 or less 
fhan 3 minutes, which means the transfer 
linction was obtained with minimum ex- 


COMPUTER CIRCUIT OF THE SYSTEM 


TRANSFER FUNCTION WITH 


STABILITY COMPENSATION ADDED. 


penditure of engineering time. The 
results of the computer program show 
that the analytic expression of the data 
is 

KG 


bes (s/0.4+1)(s/1.25+1)(841.8) 
s(s/0.25+ 1) [(s/4.3)?+(2)(0.45)s/4.34+ 1] 


For the first run, the computer program 
was set up—because of the author’s 


faulty guess—for an expression with the 
form 


a K(s/tu+1) 
3(s/ta+ i) [(s/ws1)?+2 pars/os1+ 1] 


A short study of the first set of data from 
the computer program, as a result of this 
approximation, indicated that the actual 
expression had an additional factor in the 
numerator, and this was added for the 
subsequent trials. The third and fourth 
computations which were made were 
necessary only because of lack of experi- 
ence with this program since they merely 
decreased values of the increments and of 
the allowable error. However, by using 
larger increments on the first few runs, 
the computer operating time was un- 
doubtedly shortened at the expense of 
using slightly more engineering time. 
This consideration is one that probably 
has to be worked out at each individual 
facility, being a function of available 
computer time and cost versus engineer- 
ing time and cost. 

After the analytic transfer function of 
the inherent section of the system was 
calculted, a simple lead-lag network was 
added to enable the system to operate in 
a stable mode at the open-loop gain of 
about 60 decibels. Figs. 1 and 3, respec- 
tively, show the open-loop Bode plots 
without and with the stabilizing com- 
pensating network. The transfer func- 
tion, then, of the stability-compensated 
open-loop system becomes 


KG 
_ (841.8)(s/0.4+1)(s/1.25+1)(s/31.24+1) 
~-$[(s/4.3)2+(2)(0.45)s/4.3-+1] X 
(s/0.25+1)(s/156+1) 


Fig. 5 (left). Block 
diagram of analog 
computer setup 
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VARIABLE) 


Fig. 6 (right). Tran- 
sient response of 
optimized system 
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and the closed-loop response can be cal- 


; culated to be 


KG 
(s/0. 4+-1)(s/1.25+1)(s/31.2+1) 


oF 00000165s5-+ 0.000264s4-++0.0652s8+- 
2.111s?+3.333s+1 


The inverse Laplace transform of the 
system, incorporating step and ramp in- 
puts, were obtained by appropriately 
modifying the expression for each type of 
input and solving for the output as a 
function of time. 


Checking System Response 


The calculated and empirical response 
of the system were checked and found to 
be in close agreement, but the response, 
shown in Fig. 4, was clearly not optimum. 
At this point, it must be remembered 
that these data and tests are for small 
signals which keep the system in a linear 
domain. To decrease the settling time 
of this positioning system, it seems 
obvious that a signal which would de- 
crease or eliminate oscillations and over- 
shoot is necessary. 

Now the performance-optimizing that 
was done was only partial and, as is 
usually the case, it was a compromise 
between cost, reliability, and perform- 
ance. In this case, it consisted of modi- 
fying the armature voltage with a voltage 
proportional to current, and using this 
asadamping signal. This particular part 
of the design was aided by an analog 
computer simulation of the computed 
system equation. 

The integral of the absolute error of the 
system was used as the criterion for op- 
timizing the response time. The damp- 
ing signal mentioned previously was sent 
through a variety of frequency-sensitive 
networks with variable gains and a 
minimum value of the integrated ab- 
solute error was chosen as being this 
system’s optimum performance. A sche- 


X 
7 
OUT (t) 


8 9 


wo ll 2 
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matic representation of how this was 
accomplished is shown in Fig. 5. 

The results of this study showed that, 
if the damping signal described was sent 
through a network with a transfer func- 
tion of 0.045/(0.1S+1), the integral of 
the absolute value of the error was a 
minimum. The error could be made even 
smaller by changing the system gain by 
a factor of 10 and by modifying the com- 
pensation still further, but the perform- 
ance improvement would be slight and the 
additional cost required to effect this 
change would not be warranted in this 
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case. The time response of the opti- 
mized system is shown in Fig. 6. 

The optimizing of a system, using an 
accurate mathematical representation of 
the system on a computer, as compared 
to adjusting compensating networks to 
improve performance in the field, is a 
debatable question, probably resolved 
best by doing both. 


Conclusions 


In summary, customer-engineered sys- - 


tems can be built more accurately, eco- 


nomically, and efficiently if the gen 
type of problem can be set up for 
tion on modern computers. For 
initial problem, additional time is spe 
setting up the computer program to han 
a general problem of its type, but this 
compensated for in customer satisfa: 
on the first job and the realization 
future problems in the same category wil 
be cut to a minimum. _” 
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PHE PERFORMANCE required of 
M new cranes and hoists is becoming 
gularly more demanding; specifications 
pw call for capabilities not previously 
ailable. Safety and dependability have 
ways been of primary importance in 
bisting work, and to these requirements 
e now often added a wide speed range 
¢ the rapid movement of light loads or 
he empty hook; quick, smooth accelera- 
on and deceleration of all loads with 
generative braking to reduce service on 
ie mechanical brake; effective torque 
mitation, both steady-state and tran- 
ent, to limit the mechanical strains on 
he entire structure as well as to prevent 
kcessive overloading; controlled slow 
beeds essentially independent of load; 
ad load suspension or extremely fine inch- 
g for the positioning of sensitive loads. 
l these characteristics and more must 
2 furnished within the framework of a 
ugged, simple system requiring little or 
fo maintenance. 
| Many types of hoist drive have been 


Inployed for the satisfaction of any or all 
if these requirements. Most are well 
own and it is necessary only to keep in 
hind that some make use of a-c motors 
hth various schemes for the modifica- 
jon of the a-c motor characteristics, and 
hat d-c motors, both series-wound and 
| hunt-wound, are used in other drives. A 
haracteristic of a-c motor drives is that 


} . 


ine sizing of.the motor is a function of the 


aximum speed and maximum load that 


hust be handled. In the event that a 2- 


b-1 or 3-to-1 ratio of empty hook speed 

b rated hook load speed is desired, the 

motor name-plate horsepower that must 

e furnished will be in about the same 

roportion. 

A separate high-speed motor may be 
nished as an alternative. 


The d-c 
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series motor has desirable characteristics 
for general hoisting service, but has limi- 
tations in that large currents must be 
switched or dissipated in heavy resistors, 
that operation at very light loads is un- 
certain, and that low-speed operation or 
inching at light hook loads is difficult. 
The d-c shunt motor is normally a device 
of relatively constant speed. When used 
in combination with an adjustable voltage 
supply it can provide a wide speed range 
with speed regulation nearly independent 
of load. 


The Circuit 
The hoist drive described in this paper 


makes use of a shunt-wound d-c motor in 
combination with an adjustable voltage 


300 POINT 5 


POINT 4 


HOIST 


POINT 3 


POINT 2 


PERCENT RATED 
HOOK SPEEN 


POINT 
POINT 2 


POINT 3 
100 


LOWER 


200 
POINT 4 


309 Font 5 


Fig. 1. 


10) CATA AYAVAYATAYAVAYAVAVANAVATAYAYATAVATAYAVAYAYAVAYA"AYAYAYAVAVAYAYA"AYAVAVAVAVAUAVAVAYAY) SASSO 
4 00, ) 
Eee oO RRO LGN OAR LG TAN 0 


control scheme which, when desired, im- 
parts to the motor the essential character- 
istics of the series motor, or utilizes the 
constant speed properties of the shunt 
motor under load conditions where it is 
useful. This is accomplished by the use 
of the constant-horsepower, or field- 
weakening, range of the motor’s character- 
istic to produce light hook speeds of two 
to four times the full-capacity speed. 
Typical load speed curves are shown in 
Fig. 1. The hyperbolic characteristics 
of the series motor are evident in the 
curves for the fourth and fifth positions of 
hoisting and lowering. This is operation in 
the motor field-weakened range. Rela- 
tively constant speed characteristics: of 
the armature-controlled range are em- 
ployed in the first three points in either 
direction for consistent pickup and ap- 
proach performance. The curves of Fig. 
1 illustrate the effect of torque (or cur- 
rent) limiting in that the hoisting speed 
falls to zero if the load exceeds some 
predetermined value. The load limit in 
the lowering direction is sufficiently 
higher than that in the hoisting direction 
to assure that the drive will lower under 
control any load that it is capable of lift- 
ing. 

The fundamentals of the required cir- 
cuitry are shown in the diagram of Fig. 
2(A). A-c power is almost universally the 
starting point for electrically operated 
hoist drives; it is converted to direct cur- 
rent in the rotating motor-generator set. 
The armature of the d-c generator is con- 
nected directly to the armature of the 
shunt hoist motor. An amplifier num- 
bered 2 is the exciter which drives the field 
of the d-c generator; amplifier 3 is the 
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Fig. 2. A—Fundamental circuit for d-c hoist 
drive. 
generation on deceleration 


exciter for the shunt field of the hoist 
motor; and amplifier 1 is the regulator. 

The basic circuit is that of an armature 
feedback adjustable-voltage drive in 
which the d-c armature voltage is forced 
to conform in magnitude and polarity to 
the low-power reference signal input to 
amplifier 1. This arrangement in itself 
would result in a constant speed shunt 
motor characteristic. Direct armature 
current negative feedback is supplied 
to amplifier 2, thereby forcing the output 
of amplifier 1 to be proportional to the 
armature current at. all loads. 

If amplifier 1 were linear without limit, 
the current from amplifier 1 to amplifier 
2 would necessarily increase in magnitude 
as the motor was loaded. As the negative 
current feedback acted to reduce the 
generator voltage, the slight reduction in 
armature voltage would force amplifier 1 
to even higher outputs. Amplifier 1, 
however, is arranged to saturate at a suit- 
able level and when the armature current 
has risen to the level which produces this 
saturation, the amplifier can no longer 
offset the effect of the negative current 
feedback to amplifier 2. Any further in- 
crease in load current will only serve to 
reduce the output voltage of the generator. 
The drive will stall at currents not much 
in excess of those at which the first 
amplifier saturates. 

The motor shunt field amplifier 3 is 
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B—With circuit to prevent excess re- | 


biased by the control circuit to supply 
full motor field current or an appro- 
priate weak field value, depending upon 
which characteristic it is desired to pro- 
duce. The output of amplifier 1 is fed 
into amplifier 3 with a polarity such as 
always to increase the motor field current. 
In those positions of the master switch 
in which full motor field current is 
called for anyway, this signal is over- 
ridden by a bias signal and has no effect 
on the operation of the drive. At those 
speed points for which the exciter ampli- 
fier 3 is biased to low output, the motor 
field current will remain at weak field 
value so long as the motor armature cur- 
rent is not sufficiently great to drive the 
output of amplifier 1 to a value larger than 
the bias signal on amplifier 3. As the 
load current tends to become larger than 
this value, the increasing output from 
amplifier 1 drives the motor shunt field to 
higher currents, causing the hoist to 
operate at lower speeds. 

The use of a special feedback winding 
on amplifier 2 carrying full armature 
current provides a fast, definite, current- 
limiting action. As a result, current 
overshoot (and hence peak torque) is 
kept down in magnitude and duration 
providing better protection to the elec- 
tric machinery and mechanical structure. 
Time lags in amplifiers 1 and 2 tend to 
develop instability in the system when 
operation is in the current limit range, but 
this is easily eliminated by a negative 
rate feedback around amplifier 1. 


Operation 


A consideration of the events following . 


a change in input to the drive will serve 
to clarify the operation of the circuit and 
to illuminate the desirable character- 
istics. Consider first the sequence of 
events following a change in master 
switch position from standstill or low- 
speed hoisting to a higher speed hoisting. 
A larger input signal to amplifier 1, un- 
balanced at first by any corresponding 
increase in armature voltage, will drive 
amplifier 2, and hence the generator field, 
to full-voltage values. As the generator 
voltage rises, the armature current will 
rise. This current will tend to limit itself 
through the feedback to amplifier 2 if it 
should exceed the desired current maxi- 
mum (usually 150% to 175% of rated 
armature current). 

If the hoist is lightly loaded, the arma- 
ture voltage will come up to the value 
called for, the armature current will be 
small, and the motor will settle at the 
required speed. If the hoist should be 
heavily loaded, the armature current will 


continue to be large, requiring a con) 
tinued large output of amplifier 1 in ordel 
to maintain the desired speed. If thd 
position to which the master switch hay 
been moved calls for weak field on thd 
drive motor, the bias to amplifier 3 wilt 
have been such to produce the necessary 
low output from that amplifier. Hows 
ever, during the acceleration period, thd 
high output of amplifier 1 will force thd 
output of amplifier 3 to remain at ful: 
value, and the acceleration will tak 
place at full motor field (and torque): 
If, when the appropriate armature volt! 
age is reached, the load current is 1 
the output of amplifier 3 will adjust ti 
the low value necessary to develop a higi 
hook speed. If the load on the hooki 
large and the armature current is high 
the continued higher output of amplifier - 
will maintain the motor field at or near iti 
full value. } 
The sequence of events is much t k 
same when the master switch is moveé 
from some hoisting position, or from stand 
still, or from a low-speed lowering poss 
tion to a higher speed lowering position 
Assume that the point called for req 
field weakening. Amplifier 1 will 
amplifier 2 and the generator field in thi 
direction to call for the required armatur: 
voltage, now of the opposite pola 
from the hoisting case. The high outpr 
of amplifier 1 keeps the exciter amplifier ° 
at full output, assuring full motor fie 
current during the acceleration. Tf thi 
load on the hook is large it will overh : 


As the motor-induced voltage reac e! 
and passes the value induced in the get 
erator, the armature current will go ti 
zero and reverse direction. The motori 
now regenerating. If the load is larg 


the regenerative current is large enot g: 
to counterbalance the load torque ani 
maintain a constant lowering speed. 


fier 2 is now in the direction to increai 
the induced voltage in the generata 
(since the regenerative current has 1 h 


current). The regulator amplifier 1, igh 
ever, is acting to maintain the se 
artnature voltage and, as that voltag 
rises, the output of amplifier 1 wi 
diminish to zero and reverse and increas 
with opposite polarity. If the load : 
heavy enough, this output’ will be s 
cient so that when it is fed into amplifi 
the motor field will be maintained at : 
level for slow and controlled lowering’ 
the heavy load. In this manner, th 


Marcu 1 


0/A4A0V .30-60 CYCLE 
A.C. LINES 


CUSTOMERS 
FUSED 


DISCONNECT 
SWIT 


stop START IM 


GENERATOR 
HOIST 


FIELD 
(1F2) 


(1F4) 


GENERATOR 
LOWER 


25 
(=) 
a 
4 


rs 
Ny 


HOIST MASTER SWITCH ISS 
CONTACT DEVELOPMENT 


4 


4 
ZY 


Xp Be 


XS 
x] 
X 


9 .C HOIST 


bd DDS 


tor field current is kept at or near full 
lue duringall accelerations and decelera- 
ms, and it becomes weak only when the 
d is shown to be small enough to war- 
t high-speed movement. The region 
operation in which lowering control is 
omplished through both the generator 
d and the motor field is an inherently 
stable region. This, however, is over- 
e by manipulation of the time con- 
nt of the motor field amplifier 3. 

The same delaying of the motor field 
plifier prevents a transient dip in motor 
id current during acceleration in 
lowering direction when the output 
amplifier 1 must pass momentarily 
ough zero. Hence, acceleration in the 
ering mode and with a heavy load 
accomplished without the large over- 
ot in speed that is characteristic of 
e types of hoist drives. 

Two factors are inherent which permit 
drive to lower safely under control 
oad larger than it can hoist. The in- 
‘ent mechanical efficiency of motor, 
ing, and drum arrangement is such 
t a larger hook load is necessary to 


ARCH 1961 


MOTOR 


IR 2R ISX 2SX 
ea el 0) (8) (ap Oy 


CURRENT LIMIT 
FEEDBACK WOG. 


COMMUTATOR FIELD 
eS a ba 


1 tH LSH | ML 
o4+—()— | 
LIRTA 
tH IL LSL | 
2 1M ie 
Ot 1} Ses! at 
IL IBR ; 
1 
Ue 1SS-C Len 
1H 
TO CUSTOMERS 
BRAKE CIRCUIT 1SS-D 138 
ISS-E IMF 
Veeeas O 14R 
1SS-F 
=| 
fe 
Hi [ee 
e|( 5) 
1SX 
(6) 
(6) 
COMMUTATOR N. C. HOIST asx 
inl FIELD GENERATOR 
hs + 
(T-30) 


f 


Fig. 3. Schematic diagram for hoist drive 


produce the same regenerative current as 
is produced in hoisting a given load. Be- 
cause the output of amplifier 1 must re- 
verse and build up to saturation from a 
certain “‘negative” output when lowering 
an overhauling load, more armature 
current will be required to saturate the 
amplifier. If required, the circuit can be 
easily modified to provide a further built- 
in difference in limits for each direction. 

In cases where a wide field range is 
employed (three or four to one it is possi- 
ble that during deceleration from a high 
speed in the field range, too large a volt- 
age may be built up in the motor armature 
so that saturation of the generator pre- 
vents the effective operation of the current 
limit. In this case a spillover circuit can 
be employed as shown in Fig. 2(B) to hold 
the motor field at a weak value until the 
output of amplifier 2 begins to drop off. 
This will provide a deceleration at con- 
stant armature current until the motor 
reaches about base speed; then the field 
will be permitted to increase as the motor 
comes to standstill. 

A simplified schematic of the control to 
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th 


provide this performance is presented in 
Fig. 3. Here, magnetic amplifiers are 
shown for each of the amplifiers described 
above. Rotating, electronic, or semi- 
conductor amplifiers could be used in their 
stead. The power circuit contains no 
contactors or reversing switches. The 
heavy current feedback keeps the cir- 
culating currents during standstill to low 
values. The push-pull regulator circuit 
provides a safe, gentle lowering of even 
the heaviest load if the mechanical brake 
should fail with the master switch in the 
OFF position. The reference is calling for 
zero armature voltage under such a condi- 
tion, and the lowering speed will be less 
than that associated with first-point 
lowering. 

The circuit lends itself easily to the 
use df stepless control. By the use of a 
potentiometer input and by limiting the 
reference voltage to that used for first- 
position hoisting or lowering, a contin- 
uously variable speed band between about 
10% of base speed in either direction can 
be achieved. Of course, it is then possi- 
ble to hold the load stationary on motor 


39 


K | % =0.30 OHM 
. 8&5 AY 2 k= 10 R=0.3 
— Ae Re 1 R6=100 
>) eat ea HOTOR : Tg=0.2 SEC 
Ss é ke=2.0 ——— TS pen: © 068 RAD/SEC = BASE SPEED 
Tg ke kt 6 & RAD/SEC x FIELD AMPERE 
3 Pe tia LB FT 
° Re ARMATURE AMPERE x FIELD AMPERE 
r i = 100 AMPERES ; 
ab WHEN LOWERING AT 2.5 x BASE SPEED, @ )=-170 RAD/SEC, i g=0.80 AMP 
= WHEN LOWERING AT 3.0 x BASE-SPEED, © [7-205 RAD/SEC ig =0-67 AMP 
A keA\ CIRCUIT AND CONSTANTS FOR STABILITY: CALCULATIONS FOR D.C. HOIST 
tee WiKi 
rectifier  |©2| 
: . s ele 
WAN Fig. 4. Fundamental circuit used for stability 
analysis 
2 a $$ —— 
et hex 
torque, or to inch the hook in extremely +4. Jdw/dt=kriig (5) 
small increments. During load suspen- ; 
kee2= Rete 1+ Tmd/dt) (6) 


sion, armature current flows through the 
stationary commutator. In general, even 
with rated load on the hook the current 
will be less than rated value. 

Momentary control is provided to 
prevent this condition from being acci- 
dentally maintained. 


Stability Analysis 


As previously mentioned, in the con- 
stant power range the drive is inherently 
unstable in the lowering mode, and special 
precautions must be taken to stabilize it. 
One method of stabilization is to slow 
down the response of the motor field cur- 
rent to change in the output signal of 
amplifier 1. How much damping of the 
motor field current is required to achieve 
stable operation? The necessary amount 
of stabilization depends on the field-weak- 
ening range. The problems involved can 
best be illustrated by analyzing a given 
drive. 

Fig. 4 shows the portion of the circuit 
which is significant in the constant power 
range. The circuit constants given in the 
figure are assumed; they are selected in 
such a way that the calculations will be 
simplified; however, they are realistic and 
could very well apply to a 25-horsepower 
hoist. For ease of calculation, it is also 
assumed that the amplifier time constants 
are negligible in comparison to the field 
time constant and the mechanical time 
constant. This assumption can be real- 
istic even with the use of magnetic ampli- 
fiers, if amplifiers with response approach- 
ing one half-cycle are applied. 

The circuit equations for both hoisting 
and lowering are: 


€2 = ke( e€1—e4) (1) 
e(1-+T, d/dt) =kea(e2—ki) (2) 
4 =es5 = Ri (3) 
5 = ketvig (4) 
40 
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Because products of variables are con- 
tained in the equations the equations are 
they can be linearized by 
selecting certain operating points and 
investigating the stability in the vicinity of 
Let one operating point for 
load torque be 7 and the corresponding 
Apply a 
small variable change Ar; to the load, 
which will then assume the value 77= 7+ 
The corresponding values of the 


nonlinear; 


these points. 


other variables be w, é21, @a1, etc. 


ATi. 
other variables will be w=a;+Aw, @= 
éx+ Aes, etc. 
1 through 6 yields 


Aéo = koAe; — koAea 

Aes. 1+ Td /dt) = kesAe2— 3st 

Aes— Aes = RAt 

Aes =e Ake+keinAw+keAighw 
Ari+JdAw/dt =k,i:Aigt+krig Ai+k, Aihis 
keeAeo= AigRe(1+ Td /dt) 


Working with small variations only, 
the product terms of two incremental 
variables can be dropped so that linear 
Applying the Laplace 
transform to the equations after dropping 
the product terms, and letting capital 


equations result. 


letters represent the transforms, yields 
AE2=k,AE,—koAF, 

AEA 1+T ys) = kgsAEo—kgsAl 
AFE,=AE;= RAI 

AEs = Rew Al 6+ Reig: Aw 

Ari +-JAws =k, Alg+kzig AI 
ksAEo=Re(1+ Tms)ATs 


From these equations, the response to 
an incremental change in any of the 
variable A-terms can be calculated. It 
is elected to investigate the stability by 
calculating the response in speed to an 
incremental step change in load torque 
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Substituting in equations 


Ary, assuming constant reference volt; 
Ae:=0. Substituting AZi:=0 and @ 
the constants from Fig. 4, and leavs 
the motor time constant T;, in letter for 
yields the following: 


For lowering at 21/2 times base speed 


(w;=—170; t1=100; i =0.8) 
3.0—(0.402Tm-+0.0006)s— 
Ae 0.0006 Tims? 
Ar, 3.08-++(1.37I'm—1.81)s-+ (0.4057 me 
0.0006 )s?-++0.0006 Tms* 


For lowering at 3 times base speed 


(a1 = —205; y= 100; 161 =0)53) 
3.6 —(0.402Tm-+0.0006)s — 
Ao _ 0.0006 Tins? 
Ati ~ 3.8-+(1.34Tm—2.6)s+(0.28T m4 
0.00042)s?-++-0.00042 T;,s% 


Clearly, the response of the hoist 
each of the lowering cases shown, eqt 
tions 7 and 8, will be unstable unless #1 
motor field time constant is greater th: 
a certain amount. For a stable syste 
the roots of the denominator of th 
response equations must all have neg 
tive real parts. The value of Tm 1 
quired to assure stability is found to be 
least 1.82 seconds in the first case, 
1.95 seconds in the second. 

The time constant of the motor f 
winding when operating in the field-wea 
ening range is lower than these value 
Thus, it is necessary to control the tin 
constant of the motor field exciter, amp 
fier 3. The values calculated are bor: 
out by the results of tests. As the f - 
control or constant horsepower rangi : 
extended the field circuit must be m 
to respond more slowly in order to m ai 
tain stability. This presents no 
handicap when it is working at thre ° 
four times base speed. To keep the de 
to a minimum, the higher lowering spe | 
can be achieved with the use of hig 
armature voltages necessitating a 1o} w 
field range. Since no stability proble 
exists when hoisting in the constant po 
range, the field circuit damping ca: ; | 
eliminated in this mode; when lowe r 
a variable or nonlinear damping may 
suitable to speed the over-all response : 


Summary 


Careful manipulation of control fur 
tions has made it possible to gene ra 
desirable motor characteristics utili : 
a basically standard type of adjustak 
speed drive incorporating simple : : 
rugged circuitry. Series or shunt me 
characteristics can each be employed | 
that part of the work which it fits best. 
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